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Executive  Summary 


Wetlands  are  increasingly  at  risk  from  human 
alteration  of  the  landscape.  Although  site-spe- 
cific activities  like  dredging,  ditching,  filling, 
draining,  and  plowing  have  the  most  direct  and 
obvious  impacts  on  wetland  integrity,  activities 
within  the  surrounding  catchment  can  also  lead 
to  degradation  by  changing  wetland  hydrologic 
function,  increasing  nutrient  and  sediment 
loads,  and  providing  a  conduit  for  the  spread  of 
invasive  and  exotic  species.  With  the  wide- 
spread adoption  of  GIS  technology,  it  has  be- 
come possible  to  characterize  large  landscapes 
and  identify  potential  stressors  from  existing 
datasets.  Because  so  much  information  is 
available  on  a  desktop  computer,  the  U.S.  En- 
vironmental Protection  Agency  advocates  the 
use  of  GIS-based  landscape  analysis  to  provide 
a  preliminary  assessment  of  wetland  condition 
in  a  project  area  (Level  I),  before  conduct- 
ing field-based  rapid  (Level  II)  and  intensive 
(Level  III)  assessments. 

Although  most  Level  I  assessment  approaches 
are  developed  with  best  professional  judgment, 
when  field  data  is  available,  it  can  support 
development,  calibration  and  validation  of 
metrics.  In  Montana,  while  we  do  not  have  a 
wetland  condition  dataset  providing  detailed, 
specific  measures  of  hydrologic,  physical,  and 
biotic  parameters,  we  do  have  rapid  assessment 
data  on  over  a  thousand  wetlands  across  the 
state.  Our  goal  in  this  study  was  to  determine 
whether  we  could  use  this  data  to  identify 
landscape-level  metrics  with  a  good  ability  to 
predict  wetland  condition,  or,  at  the  least,  to 
calibrate  and  validate  a  best  professional  judg- 
ment-based tool. 

From  a  review  of  the  literature,  we  identified 
a  number  of  landscape-scale  metrics  that  are 
widely  believed  to  influence  wetland  condi- 
tion. We  calculated  values  for  these  metrics  in 
several  different  buffer  distances  for  a  random 


sample  of  591  wetlands.  We  ran  an  ANOVA 
to  determine  the  metrics  that  had  significant 
relationships  to  the  field-determined  overall 
condition  scores.  Each  significant  metric  was 
further  examined  through  a  linear  regression 
to  determine  the  strength  of  the  relationship. 
Metrics  were  also  combined  and  analyzed 
in  multiple  stepwise  regression  routines  to 
evaluate  interactions,  and  assessed  with  CART 
(Classification  and  Regression  Tree)  analysis. 

At  the  6th  code  Hydrologic  unit  (HUC),  1 
kilometer,  500  meter,  and  200  meter  buffer  dis- 
tance, the  combined  metrics  of  percent  forest 
cover,  road  density,  and  number  of  stream  road 
crossings  had  the  strongest  predictive  value  for 
overall  score.  However,  the  R-squared  value 
for  the  three  metrics  combined  was  only  0.14 
at  the  6th  code  HUC  level.  When  measured 
at  the  buffer  level,  the  combination  of  metrics 
had  an  even  lower  R-squared  value.  At  1 ,000 
meters,  the  R-squared  value  dropped  to  0.13, 
at  500  meters  it  was  0.11,  and  at  200  meters, 
the  R-squared  value  was  0.09,  and  road  density 
was  no  longer  significantly  correlated  with 
overall  score.  We  had  observed  that  there  was 
a  strong  ecoregional  skew  in  the  condition 
scores,  with  wetlands  in  mountain  ecoregions 
having  a  higher  average  score  than  wetlands 
in  plains  ecoregions  (0.93  vs.  0.77  on  a  scale 
of  0.0  to  1 .0).  Therefore,  we  split  the  assess- 
ment data  into  a  mountain  and  a  plains  subsets 
and  reran  the  analysis.  With  the  data  divided, 
percent  forest  was  no  longer  significant  at  any 
scale.  For  wetlands  in  the  mountain  ecore- 
gions (n=262),  road  density  was  the  only  met- 
ric that  was  significant  at  all  levels,  although 
the  R-squared  value  was  never  higher  than 
0.07.  In  the  1  kilometer  buffer,  the  percent- 
age of  crop  agriculture  was  also  significant, 
although  it  had  no  significance  at  other  buffer 
distances.  In  the  plains  ecoregions,  no  metrics 
were  significant  at  200  meters.  Percent  natural 
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grassland  and  road  density  within  500  meters 
were  both  significantly  correlated  with  overall 
score  but  had  very  low  R-squared  value  (0.02 
and  0.01,  respectively).  At  the  1,000  meter 
buffer  scale,  only  the  number  of  stream  road 
crossings  was  significant.  No  metric  was  sig- 
nificantly correlated  to  overall  wetland  condi- 
tion when  measured  at  the  6th  code  HUC  level 
in  either  the  mountain  dataset  or  the  plains  da- 
taset.  When  we  added  an  environmental  vari- 
able (relative  effective  armual  precipitation)  to 
the  analysis,  we  found  it  had  high  predictive 
value  for  the  dataset  as  a  whole,  and  within  the 
subset  of  mountain  ecoregions.  In  the  plains, 
where  it  varied  less,  it  was  not  significant. 

Using  best  professional  judgment,  we  then 
built  a  Montana  Landscape  Integrity  Model 
(MT-LIM)  and  used  the  dataset  to  calibrate  it. 
An  additional  1 00  points  drawn  from  the  initial 
dataset  and  1 80  points  assessed  specifically 
for  this  study  were  used  in  validation.  The 
model  is  an  inverse  weighted  distance  model 
premised  on  the  idea  that  ecosystem  processes 
and  functions  achieve  their  fullest  expression 
in  areas  where  human  activities  have  the  least 
impact.  In  the  case  of  wetlands,  it  presumes 
that  wetland  condition  will  be  poorest  in  close 
proximity  to  roads,  commercial  or  industrial 
development,  urban  areas,  resource  extraction 
sites,  or  hydrologic  modifications.  The  model 
is  built  into  a  single  raster  layer  covering  the 
entire  state.  The  raster  pixel  size  is  30  meters 
by  30  meters,  or  900  square  meters.  Pixel  val- 
ues range  from  100  to  745. 

The  model  was  used  to  calculate  a  mean 
landscape  integrity  score  for  pixels  within  100 
meters  of  a  wetland.  This  score  was  combined 


with  a  relative  effective  precipitation  value 
from  the  assessment  point,  and  wetlands  were 
assigned  to  an  ordinal  condition  class  (A,  B,  C 
or  D)  using  thresholds  we  identified  through 
calibration.  When  compared  to  the  condition 
classes  assessed  in  the  field  using  rapid  as- 
sessment methods,  this  approach  accurately 
predicted  the  measured  condition  in  50-55%  of 
cases,  depending  on  field  method.  In  mountain 
ecoregions,  it  accurately  predicted  A-ranked 
wetlands  in  75%  of  the  cases.  Sixty-five  per- 
cent of  the  plains  wetlands  and  83.5%  of  the 
mountain  wetlands  were  classified  acceptably, 
i.e.,  the  classifications  were  either  accurate  or 
no  more  than  one  rank  higher  than  what  was 
assigned  in  the  field. 

This  study  demonstrated  the  potential  of 
landscape-level  metrics  and  models  to  pre- 
dict wetland  condition  using  remotely-sensed 
data  in  Montana.  At  the  same  time,  it  showed 
that  environmental  variables  and  site-specific 
activities  may  be  far  more  important  drivers 
of  wetland  condition  than  land  uses  occurring 
at  a  broader  scale.  Human  disturbances  can 
affect  wetlands  at  multiple  scales,  and  because 
they  tend  to  interact  both  with  each  other  and 
with  environmental  variables,  it  is  exceedingly 
difficult  to  quantify  or  predict  their  effects.  Ul- 
timately, landscape-level  wetland  assessment 
may  never  be  a  satisfactory  substitute  for  field 
assessments.    Nonetheless,  when  field  assess- 
ments are  not  practical,  landscape-level  as- 
sessment is  the  only  alternative.  With  comple- 
tion of  the  Landscape  Integrity  Model  and  an 
approach  for  integrating  it  with  precipitation 
data,  we  now  have  a  tool  offering  a  quick  and 
efficient  way  to  begin  that  task. 
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Introduction 


Wetlands  are  increasingly  at  risk  from  human 
alteration  of  the  landscape  (Allan  2004,  Hychka 
et  al.  2007).  While  site-specific  activities  like 
dredging,  ditching,  filling,  draining,  and  plowing 
are  the  most  obvious  examples,  alterations  that 
occur  within  the  surrounding  catchment  can 
negatively  impact  wetlands  by  changing  wetland 
h\drologic  function  (Gregory  et  al.  199! ), 
increasing  nutrient  and  sediment  loads  (Skagen  et 
al.  2008),  and  providing  a  conduit  for  the  spread 
of  invasive  and  exotic  species  (Magee  and  Kentula 
2005).  Consequently,  most  wetland  assessment 
protocols  incorporate  metrics  that  attempt  to 
capture  the  effects  of  landscape-level  impacts, 
ty  picallv  by  measuring  the  extent  and  condition  of 
the  surrounding  buffer  or  by  including  a  checklist 
of  stressors  occurring  in  the  nearby  landscape 
(Faber-Langendoen  et  al.  2008). 

With  the  widespread  adoption  of  GIS  technology, 
it  has  become  possible  to  characterize  large 
landscapes  and  identify  potential  stressors  from 
existing  datasets  (LaGory  and  Kuiper  2004, 
Phillips  et  al.  2005).  The  USGS  National  Land 
Cover  Dataset  (NLCD),  derived  from  30-meter 
Landsat  TM  imagery,  classifies  land  cover  and  land 
use  into  16  categories  (Table  I ),  allowing  rapid 
evaluation  of  the  proximity  of  target  wetlands  to 
known  disturbances  like  agriculture  and  urban 
development.  By  combining  land  cover  data  with 
similar  datasets  that  identify  and  depict  roads, 
water  features,  and  topography,  researchers  have 
produced  synoptic  maps  that  broadly  predict 
wetland  impairment  (Brooks  et  al.  2004).  This 
success  has  contributed  to  the  growing  popularity 
of  a  three-tiered  approach  to  wetland  assessment 
(Kentula  2007).  In  this  approach,  readily 
available  GIS  datasets  are  used  to  create  an  initial 
characterization  of  wetland  condition  based  on 
landscape-level  indicators  (Level  I).  A  subset  of 
wetlands  assessed  at  Level  I  can  then  be  fiirther 
evaluated  in  the  field  using  rapid  assessment 
protocols  (Level  II),  and  a  subset  of  those  can  be 
subjected  to  intensive  assessments  using  metrics 
of  floristic  quality,  soil  status,  hydrology,  and/or 
biotic  richness  (Level  III).  If  Level  III  metrics 
and  scoring  thresholds  have  been  derived  from 


I'ahle  I.  National  Land  Cover  Datasel  (2001)  cover  classes. 


CODE 

CLASS 

11 

Open  water 

12 

Perennial  ice  or  snow 

21 

Developed,  open  space 

22 

Developed,  low  intensity 

23 

Developed,  medium  intensity 

24 

Developed,  high  intensity 

31 

Barren  land 

41 

Deciduous  forest 

42 

Evergreen  forest 

43 

Mixed  forest 

52 

Scrub/shrub 

71 

Grassland/herbaceous 

81 

Pasture/Hay 

82 

Cultivated  crops 

90 

Woody  wetlands 

95 

Emergent  herbaceous  wetlands 

measured  parameters  in  wetlands  of  known 
condition,  then  scores  from  a  Level  III  assessment 
can  be  used  to  validate  and  calibrate  the  other 
assessment  levels  (Hychka  et  al.  2007).  In  the 
absence  of  these  scores  or  data,  a  more  general 
Level  I  tool  can  be  crafted  using  indicators  derived 
from  best  professional  judgment  (Brooks  et  al. 
2004)  or  from  indicators  that  have  been  used 
successfijlly  in  other  states  (Mack  2007). 

However,  existing  Level  I  approaches  may  not 
be  easily  applied  to  Montana.  For  the  most 
part,  these  tools  were  initially  developed  in  the 
eastern  U.S..  where  natural  landscapes  are  more 
homogenous.  In  Pennsylvania,  for  example, 
percent  natural  forest  within  a  wetland  buffer  is 
a  good  indicator  of  wetland  disturbance,  because 
the  presettlement  norm  was  continuous  forest;  the 
absence  of  forest  is.  by  definition,  a  sign  that  the 
landscape  has  been  altered  (Brooks  et  al.  2004, 
Hychka  et  al.  2007).  In  Montana,  by  contrast,  there 
is  no  single  presettlement  vegetation  class,  and 
in  the  case  of  some  classes  such  as  grassland  and 
shrubland,  human  disturbances  typically  manifest 
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as  changes  in  condition  and  quality,  rather  than  in 
fragmentation  and  patchiness.  Metrics  based  on 
agriculture  present  similar  problems.  In  the  Great 
Lakes  region,  where  agricultural  development  has 
been  found  to  have  measurable  effects  on  wetland 
condition  (Brazner  et  al.  2007),  agriculture  is 
mostly  crop-based,  and  is  discernible  even  by 
the  naked  eye  on  satellite  imagery.    In  western 
rangelands,  where  most  of  our  agricultural  areas  are 
devoted  to  free-ranging  livestock  rather  than  crops, 
it  is  much  more  difficult  to  pinpoint  concentrated 
impacts  without  intensive  image  analysis  and 
post-processing.  Neither  the  NLCD  nor  the  more 
detailed  ReGAP  vegetation  maps  distinguish 
rangeland  land  use  from  herbaceous  land  cover. 
Consequently,  areas  with  heavy  livestock  use  and 
areas  where  livestock  are  absent  or  excluded  are 
both  classified  as  "grassland/herbaceous,"  generally 
considered  to  be  a  "natural"  land  cover  class  in 
Level  I  assessments.  Moreover,  even  if  we  could 
identify  rangelands  with  heavy  current  grazing 
pressure,  we  would  not  necessarily  identify  all 
areas  affected  by  grazing.  Many  landscape-level 
grazing  impacts,  such  as  loss  of  riparian  vegetation 
and  downcutting  or  entrenchment  of  streams,  have 
occurred  in  the  past  and  persist  even  in  the  absence 
of  livestock.    Finally,  in  rural  states  like  Montana, 
with  its  large  land  area  and  small  population 
footprint,  the  assumptions  of  Level  I  assessments 
may  not  apply.  For  example,  road  proximity  or 
road  density  is  typically  assumed  to  be  a  stressor 
on  wetlands.  However,  the  effect  of  these  presumed 
stressors  is  probably  very  different  in  a  state 
where  over  a  hundred  thousand  vehicles  travel  a 
highway  each  day  and  in  a  state  where  a  highway 
may  support  less  than  a  thousand  vehicles  a  day. 
Because  of  lower  levels  of  intensity,  landscape- 
level  factors  may  not  have  the  same  impacts  on 
wetland  condition  in  rural  states  as  in  urban  ones. 


and  may  generally  be  less  significant  than  direct, 
site-specific  factors  like  plowing,  excavation,  and 
damming,  or  natural  stressors  like  drought. 

Nevertheless,  Level  I  assessment  tools  hold  great 
appeal  in  Montana,  where  the  size  and  diversity 
of  the  landscape  are  such  that  even  probabilistic, 
rapid  assessments  will  take  years  to  execute  on  a 
statewide  basis.  If  the  assumption  that  landscape- 
level  stressors  compromise  wetland  integrity  is 
correct,  then  Level  I  assessments  could  be  used 
to  identify  areas  where  these  stressors  are  at  play, 
thereby  helping  resource  managers  identify  areas 
where  more  detailed  monitoring  and  assessment 
projects  should  be  done.  Similarly,  a  Level  1 
assessment  tool  could  help  delineate  areas  with  few 
threats  and  identify  where  high  quality  wetlands, 
and  potentially  uncommon  species,  are  likely  to  be 
found. 

Therefore,  we  wanted  to  determine  whether  it 
was  possible  to  develop  a  Level  I  assessment  tool 
for  Montana  and  to  base  it  on  a  combination  of 
best  professional  judgment,  metrics  derived  from 
the  literature,  and  field  data.  Although  Montana 
is  in  the  early  stages  of  conducting  Level  III 
assessments  across  the  state,  the  Montana  Natural 
Heritage  Program  has  completed  over  a  thousand 
rapid  assessments  of  wetland  condition  in  the 
course  of  amphibian  surveys  using  a  protocol 
developed  by  the  Montana  Department  of 
Environmental  Quality  (MTDEQ  2005).  We  chose 
these  data  as  our  starting  point  for  determining 
which  identifiable  landscape-level  stressors  might 
be  correlated  with  a  loss  of  wetland  condition.  If 
we  could  identify  these  stressors,  then  it  would 
be  possible  to  build  a  Level  I  assessment  tool 
appropriate  for  Montana  and  provide  a  template  for 
other  rural  western  states  to  adapt. 


Methods 


The  Dataset 

The  wetland  assessment  data  used  in  the  initial 
analysis  and  model  calibration  were  collected  by 
Montana  Natural  Heritage  Program  amphibian 
survey  crews  between  2004  and  2007.  The 
amphibian  survey  sampling  scheme  delineates 
1 !  geographic  strata  based  on  ecoregions  and 
subbasins.  Within  those  1 1  strata,  there  are  up 
to  three  land  ownership  strata  (>40%  public, 
>40%  tribal.  <40%  public  or  tribal).  In  total, 
there  are  28  ""target  populations"  from  which 
subvvatersheds  (USGS  6th  code  hydrologic  units) 
can  be  randomly  selected  in  numbers  proportionate 
to  the  total  area  and  number  of  subwatersheds 
in  the  sampling  frame.  Everv'  lentic  wetland  in 
the  selected  subwatershed  is  then  visited  and 
assessed  using  a  protocol  developed  by  the 


Montana  Department  of  Environmental  Quality 
(MTDEQ).  This  study  design,  while  appropriate 
for  amphibian  breeding  site  occupancy  surveys, 
was  not  entirely  appropriate  for  our  analysis,  as 
many  of  the  wetlands  in  any  given  subwatershed 
were  in  very  close  proximity  to  each  other  (<100m) 
and  were  therefore  subject  to  the  same  suite  of 
stressors.  To  minimize  autocorrelation,  we  overlaid 
a  grid  of  public  land  survey  sections  (one  mile 
by  one  mile)  on  the  sampled  subwatersheds  and 
randomly  selected  one  wetland  assessment  site 
per  square  mile  section  for  our  analysis.  This 
selection  provided  591  wetland  assessments, 
broadly  distributed  across  the  state  but  relatively 
concentrated  in  less-developed  areas  (Figure  1 ). 
We  drew  another  100  assessments  from  the  entire 
dataset  to  use  for  testing. 
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Figure  I.  Initial  datapoints  (n=59l). 


The  MTDEQ  rapid  assessment  method 
(MTDEQ-RAM)  evaluates  both  stressors  and 
condition  within  a  wetland  and  its  100  meter 
buffer.  Stressors  occurring  between  100  and  500 
meters  are  noted  but  not  scored.  Metrics  cover 
hydrogeomorphology,  vegetation,  buflFer  condition, 
and  water  quality  (Appendix  A).  Scores  on 
individual  metrics  are  combined  into  an  overall 
score  ranging  from  0.0  -1 .0.  Because  our  goal  was 
to  build  a  statewide  Level  I  tool  that  could  identify 
broad  condition  classes,  our  analysis  focused  on  the 
multi-metric  synthesis  represented  by  the  overall 
score. 

Data  Preparation 

Landscape- level  stressors  operate  at  multiple 
scales,  and  different  wetland  indicators  react 
in  different  ways.  In  a  review  of  the  literature. 
Brazner  et  al.  (2007)  observed  that  stressors  like 
row  crops  and  pastures  appear  to  operate  at  the 
watershed  scale,  while  effects  of  urbanization  are 
most  pronounced  at  the  site  level.  They  also  noted 
that  response  variables  are  not  uniformly  sensitive 
to  watershed-scale  versus  local  disturbances. 
For  example,  water  qualify  metrics  respond  to 
watershed  or  regional  disturbances  while  wetland 
vegetation  tends  to  reflect  more  localized  impacts. 

There  is  little  consistency  in  the  scale  at  which 
Level  I  assessments  are  conducted.  Mita  et  al. 
(2007)  used  a  300  meter  buffer  in  the  Prairie 
Pothole  Region  of  central  North  Dakota  because 
of  the  irregular,  hummocky  landscape.  Brooks 
et  al.  (2004)  used  a  one  kilometer  buffer,  as  do 
researchers  using  the  Landscape  Development 
Intensify  Index  developed  by  Brown  and  Vivas 
(2005;  Mack  2007).  Hychka  et  al.  (2007)  evaluated 
five  spatial  scales:  a  100  meter  buffer,  a  250 
meter  buffer,  the  upstream  catchment  area  and  the 
intersection  of  the  upstream  catchment,  and  each 
of  the  two  buffer  circles.  They  concluded  that 
the  most  predictive  metrics  were  measured  at  a 
250  meter  and  upstream  catchment  scale.  Given 
this  variabilify.  we  decided  to  calculate  metrics  at 
several  different  scales:  200  meters,  500  meters, 

'  http://gisportal.mt.gov/Portal/ 

'  Unfortunately,  one  dataset  that  we  would  have  liked  to  use,  the  National  Wetlands  Inventory,  is  not  available  in  digital  format 
for  all  of  Montana,  so  we  were  unable  to  incorporate  it  into  our  statewide  analysis. 


one  kilometer,  and  the  entire  6th  code  HUC.  To 
obtain  the  evaluation  areas,  we  buffered  each 
assessment  point  by  each  scale  using  ArcGIS  9.2 
(ESRI  2006). 

We  were  fortunate  that  Montana  is  rich  in 
most  fypes  of  GIS  data.  The  Natural  Resource 
Information  Service,  part  of  the  Montana  State 
Library,  maintains  multiple  datasets  and  makes 
them  freely  available  for  downloading  from  the 
Montana  GIS  Portal'.  Many  of  these  digital 
datasets  have  been  derived  from  national  datasets 
(e.g.,  the  National  Land  Cover  Dataset,  the 
National  Hydrography  Dataset,  the  TIGER  roads 
data,  the  National  Elevation  Dataset)  while  others 
are  Montana-specific  products.  In  our  initial 
evaluations,  we  used  only  data  that  were  broadly 
available:  land  cover,  roads,  and  hydrology". 
Although  our  primary  goal  was  to  create  a  Level 
I  tool  that  would  serve  the  needs  of  Montana  land 
managers  and  scientists,  we  were  hopeful  we  could 
find  an  approach  that  could  also  be  used  in  other 
Rocky  Mountain  states.  Therefore,  in  the  initial 
analysis  we  did  not  include  data  that  might  not 
exist  in  all  states,  such  as  detailed  cadastral  layers 
and  oil  and  gas  well  locations. 

After  buffering  the  data  points,  we  calculated 
land  cover  and  stressor  metrics  for  each  buffer 
using  the  Analytical  Tools  Interface  for  Landscape 
Assessment  3.0  (ATtlLA)  (Ebert  and 
Wade  2000),  an  EPA-developed  program  that 
runs  in  Arc  View  3.3  (ESRI  2000)  with  the  Spatial 
Analyst  extension.  Results  were  exported  to 
ArcGIS  9.2  for  further  analysis  and  data  display. 
We  ran  an  ANOVA  to  determine  the  metrics  that 
had  significant  relationships  to  the  overall  scores 
at  each  buffer  distance.  Each  significant  metric 
was  further  examined  through  a  linear  regression 
to  determine  the  strength  of  the  relationship. 
Metrics  were  also  combined  and  analyzed  in 
multiple  stepwise  regression  routines  to  evaluate 
interactions.  All  analyses  were  performed  in 
SYSTAT  1 1 .2.  (Systat,  Inc.)  Table  2  shows  the 
calculated  metrics  used  in  this  stage  of  the  analysis. 


Tahlt'  2    Metrics  used  in  initial  analvsis. 


METRICS  BY 

UNIT  OF 

CATEGORY 

MEASUREMENT 

Land  cover  and  land  use 

Percent  forest 

percent 

Percent  grassland 

percent 

Percent  shrubland 

percent 

Percent  total  agriculture 

percent 

Percent  crop  agriculture 

percent 

Percent  pasture/hay 

percent 

Percent  urban 

percent 

Percent  natural 

percent 

Road  Densitv 

km/km - 

Patch  density  and  size  metrics 

Mean  forest  patch  Size 

hectares 

Mean  grassland  patch  size 

hectares 

Number  of  forest  patches 

count 

Number  of  grassland  patches 

count 

Forest  patch  densit> 

count/  km- 

Grassland  patch  density 

count/  km- 

Hydrologj  metrics 

Percent  altered  streams  in  buffer 

percent 

Stream  road  crossings 

count 

Reservoir  densit>' 

count/  km- 

Water  rights  points  of  use 

count/  km- 

Population  and  ownership  metrics 

Population  density 

count/  km- 

Percent  public  land  ownership 

percent 

Percent  land  in  managed  status  (wilderness. 

RNA.  conservation  easements)           percent 

Percent  land  owned  by  timber  companies 

percent 

When  the  initial  analysis  failed  to  yield  sufficient 
predictors  for  model  development,  we  added  an 
environmental  variable,  relative  effective  annual 
precipitation  (REAP),  to  the  analysis  (NRCS 
2008).  REAP  differs  from  standard  annual 
precipitation  measurements  in  that  it  represents 
the  amount  of  moisture  available  at  a  location, 
accounting  for  precipitation,  slope,  aspect,  and 
soil  properties.  We  used  a  simple  least  squares 
regression  of  overall  score  against  the  interpolated 
value  for  REAP  at  the  wetland  assessment  point, 
using  the  entire  dataset.    This  variable  had  good 
predictive  value  for  modeling  wetland  condition,  so 
we  incorporated  REAP  into  subsequent  analyses. 
We  then  used  CART  (Classification  and  Regression 
Tree)  analysis  to  further  analyze  data  for  model 
development.  CART  is  a  nonparametric  method 


frequently  used  to  explore  ecological  response 
data  involving  multiple  predictor  variables  that 
may  be  interacting  in  ways  that  are  difficult  to 
model  with  linear  approaches  (Brazner  et  al.  2007, 
De'Ath  and  Fabricus  2000).  Classification  trees 
help  to  pinpoint  predictor  variables  that  maximize 
within-group  homogeneity  in  a  dataset.    The 
CART  analysis  was  performed  using  SYSTAT 
1 1 .2"s  TREE  program,  which  uses  an  automatic 
interaction  detection  algorithm  to  find  the  best 
fit  for  the  regression.  The  initial  cluster  of 
values  is  subjected  to  stepwise  splitting  using 
all  possible  predictor  variables,  seeking  the 
variable  values  that  minimize  within-group  sums 
of  squares.  Each  group  is  repeatedly  split  until 
the  proportional  reduction  in  error  (PRE)  falls 
below  a  set  percentage  based  on  a  least  squares 


loss  function.  CART  works  best  when  there  is  a 
categorical  response  variable,  so  we  changed  the 
numeric  overall  scores  to  letter  ranks,  setting  A  at 
0.9  or  above,  B  at  0.75-0.89.  C  at  0.50-0.74,  and  D 
at  0.49  or  below'.    CART  represents  interactions 
between  predictor  variables  as  branches  from  a 
split  point,  and  continues  to  split  data  along  each 
branch  until  a  predefined  goodness-of-fit  criteria 
is  not  met,  and  the  branch  concludes  in  a  terminal 
node.  The  process  ends  when  all  branches  have 
terminal  nodes,  at  which  point  users  can  "prune" 
the  tree  back  to  a  specified  number  of  terminal 
nodes.  CART  analyses  were  run  on  data  at  three 
different  buffer  distances:  200  meters,  500  meters, 
and  1 ,000  meters. 

Model  Development 

Our  intent  was  to  use  the  results  of  the  CART 
analysis  to  select  metrics  for  a  model  that 
would  capture  condition  classes  across  the  state. 
However,  the  analysis  did  not  yield  enough  metrics 
to  do  this.  At  this  point,  there  were  three  major 
challenges: 

1 )  The  existing  dataset  did  not  completely 
capture  the  range  of  landscape-level  stress- 
ors that  exist  across  the  state.  For  example, 
while  there  was  a  statistically  significant 
relationship  between  both  crop  agriculture 
and  urbanization  and  wetland  condition, 
those  relationships  explained  less  than  1%  of 
the  variability  in  the  data,  simply  because  so 
few  wetlands  in  the  dataset  were  near  urban 
or  agricultural  areas. 

2)  Similarly,  the  stressors  affecting  some  of 
the  low-scoring  wetlands  were  not  identifi- 
able using  the  NLCD  alone.  This  was  in 
part  due  to  the  age  of  the  imagery  used  in 
the  NLCD;  several  areas  that  are  classified 
as  forested  have  been  intensively  logged  in 
the  past  few  years.  When  we  included  land 
ownership  as  a  variable  in  our  CART  analy- 
sis, ownership  by  timber  companies  was  a 
splitting  factor. 

3)  In  all  analyses  of  our  dataset  at  every  buf- 
fer distance,  the  effect  of  environmental  vari- 
ables far  outweighed  the  impacts  of  human 

^  The  MTDEQ  RAM  form  assigns  categorical  ranks  at  different  cutoff  points.  We  experimented  with  several  thresholds  in 
CART  before  concluding  that  these  yielded  the  clearest  results.  These  thresholds  agree  with  those  chosen  by  Faber-Langendoen 
et  al.  (2006)  in  their  ecological  integrity  assessments. 


disturbance.  REAP  appeared  as  the  driving 
factor  determining  wetland  condition. 

One  common  goal  of  wetland  assessment  and 
monitoring  programs  is  to  evaluate  the  success 
of  management  efforts  (Fennessy  et  al.  2004, 
Kentula  2007).  By  implication,  wetland  assessment 
tools  should  not  incorporate  both  environmental 
variables  such  as  precipitation  and  human 
disturbance  variables  like  land  use  because  doing 
so  makes  it  difficult  to  determine  what  management 
actions  could  improve  wetland  condition. 
However,  several  of  our  human  disturbance 
mefrics,  though  statistically  significant,  did  not 
have  strong  enough  relationships  with  wetland 
condition  to  support  a  model. 

As  noted  earlier,  most  Level  I  assessment  tools 
are  expert  systems  developed  by  application  of 
best  professional  judgment  (Tiner  2005,  Weller 
et  al.  2007).  Therefore,  we  decided  to  construct 
a  landscape  integrity  model  (LIM)  incorporating 
both  metrics  derived  from  the  analysis  and  our 
knowledge  of  landscape- level  stressors.  Landscape 
integrity  models  have  been  used  as  a  stratification 
tool  in  development  of  fine-scale  assessment 
methods  in  the  Rocky  Mountain  West  (Rocchio 
2007).  We  thought  that  such  a  model,  with  scoring 
thresholds  calibrated  to  existing  assessment  data, 
would  meet  our  primary  goals  of  identifying  areas 
where  high-quality  wetlands  were  likely  to  be 
found  and  areas  where  wetlands  might  be  most  in 
need  of  management  actions. 

Our  model  used  an  inverse-weighted  distance 
approach  (Tuflfly  and  Comer  2005)  to  create  land 
cover  and  land  use  layers  in  a  GIS.  We  used 
ANOVA  and  CART  to  evaluate  each  new  layer 
against  the  assessment  data.  CART  was  also  used 
to  set  weights  and  scoring  thresholds.  In  doing 
so,  we  attempted  to  optimize  use  of  layers  derived 
from  readily  available  datasets  (e.g.,  land  cover, 
roads,  hydrology)  while  minimizing  use  of  layers 
derived  from  datasets  that  might  not  be  available 
in  all  Rocky  Mountain  states  (e.g.,  cadastral,  mines, 
energy.  Clean  Water  Act  Section  404  permits). 


Table  3  shows  the  data  layers  included  in  model 
development.  Appendix  B  contains  a  detailed 
description  of  how  the  model  was  developed  and 
how  individual  data  layers  were  weighted  and 
scored. 

Because  REAP  was  seen  as  a  major  driver  behind 
wetland  condition,  it  cannot  be  ignored  in  Level  I 
assessments.  To  separate  environmental  impacts 
from  variables  that  are  responsive  to  management 
action,  we  devised  a  scoring  scheme  that  uses 
the  final  model  in  conjunction  with  REAP  data. 

Table  3.   G/S  data  sources  used  in  this  report. 


Model  Calibration,  Testing  and 
Validation 

We  calibrated  the  model  to  the  MTDEQ-RAM 
scores  in  the  initial  59]  datapoints  using  raw 
numeric  scores.  The  raster  layer  constituting  the 
LIM  assigns  a  score  to  each  30-meter  by  30-meter 
pixel.  To  assign  LIM  scores  to  a  wetland,  we  used 
Spatial  Analyst's  zonal  statistics  function  to  derive 
mean  scores  for  buffered  areas  around  each  point 
(100,  300  and  1.000  meters).  We  used  the  highest 
R-squared  value  to  identify  the  optimal  buffer 


TYPE 

SOURCE 

DATE 

CELL  SIZE/ 
RESOLUTION 

SOURCE  DATA 

Land  Use 

NLCD 

2001 

30  m 

LandSATTM(1993- 
2000) 

Hydrology 
(Streams,  rivers, 
lakes,  reservoirs, 
hydrologic  units) 

NHD 

2007 

1 :24,000 

1:24k  USGS 
topographic  maps 

Wetlands 

USFWS 

Varies 

1:24,000 

National  wetlands 
inventory 

Digital  Elevation 
Models 

USGS 

Varies 

30  m 

1:24k  USGS 
topographic  maps 

Roads 

TIGER/US 
Census  Bureau 

2001 

1:100,000 

1:100k  USGS 
topographic  maps 

Relative 

Effective  Annual 
Precipitation 
(REAP) 

NRCS 

2008 

10m 

DEMs,  DAY  MET 

Montana 
cadastral  data 

MT  Base  Map 
Service  Center 

2008 

N/A 

BLM  GCDB 
(boundaries) 

Montana  Public 
Stewardship 

MTNHP 

2008 

1:100,000 

Varies 

Water  Rights 
Points  of  Use 

DNRC 

Varies 

1:100,000 

BLM  GCDB 

(boundaries) 

Population 
Density 

Montana  State 
Library 

2002 

1  km 

US  Census  Bureau 

Oil  and  gas  wells 

NRIS 

Varies 

Varies 

BLM,  MBMG, 
MDNRC 

Abandoned  mines 

MTDEQ 

2007 

1/2  mile 

MTDEQ 

CWA  404  permits 

NRIS 

2003 

N/A 

ACOE 

area  within  which  LIM  scores  should  be  averaged. 
We  then  used  categorical  MTDEQ-RAM  ranks 
to  define  thresholds  for  each  LIM  rank,  adjusting 
values  until  the  highest  accuracy  (%  of  modeled 
and  measured  ranks  in  agreement)  was  achieved. 
Next,  we  tested  the  calibrated  model's  agreement 
with  ranks  derived  from  the  scores  that  had  been 
assigned  to  the  100  test  set  data  points  in  the  field. 

We  validated  the  model  on  an  independent  set  of 
wetlands  (n=l  80)  assessed  as  part  of  this  project 
during  2007  and  2008.  These  180  wetlands 
consisted  of  60  wetlands  chosen  from  a  spatially- 
distributed,  randomly  selected  set  of  1,800 
wetlands  drawn  fi-om  all  mapped  natural  palustrine 
and  lacustrine  wetlands  in  the  state.  The  remaining 
120  wetlands  were  part  of  a  targeted  sample 
compiled  to  incorporate  more  environmental  and 
disturbance  gradients  and  to  represent  areas  of 
the  state  that  are  not  mapped.  We  also  ran  the 
model  on  the  1,800  wetland  dataset  as  a  "dry  run" 
of  a  Level  I  assessment  procedure.  To  determine 
whether  the  model  could  predict  wetland  condition 


measured  with  other  Level  II  assessment  protocols, 
we  compared  it  against  a  dataset  of  1 08  Ecological 
Integrity  Assessments  (EIAs),  a  rapid  field-based 
protocol  we  are  currently  testing  in  Montana 
(Appendix  C). 

Model  accuracy  was  assessed  by  calculating 
the  percentage  of  wetlands  that  were  correctly 
classified  into  each  rank  category,  using  the  LIM 
score  as  the  predicted  score  and  the  MTDEQ-RAM 
score  or  EIA  score  as  the  measured  score.  Model 
acceptability  was  determined  by  calculating  the 
percentage  of  wetlands  that  were  either  correctly 
classified  or  classified  as  being  one  rank  higher 
than  measured  in  the  field  with  the  MTDEQ- 
RAM  or  EIA.  We  reasoned  that  site-specific 
factors  affecting  the  wetland  could  account  for  the 
difference  if  the  measured  condition  was  one  rank 
lower  than  the  predicted  condition.  However,  if 
measured  condition  was  one  or  more  ranks  higher 
or  two  ranks  lower  than  predicted  condition,  this 
was  considered  to  be  unacceptable  error. 


Results 


Initial  Analysis:  Step  1 

The  591  wetland  data  points  used  in  the  initial 
analysis  were  spread  across  five  ecoregions 
(Omemik  1987):  the  Northwestern  Great  Plains 
(n=207),  the  Northwestern  Glaciated  Plains 
(n=122),  the  Northern  Rockies  (n=125),  the 
Middle  Rockies  (n=108).  and  the  Canadian 
Rockies  (n=29).  CART  anaK  sis  showed  a  distinct 
ecoregional  break  in  the  data,  w  ith  the  three 
mountain  ecoregions  having  significantly  higher 
scores  than  the  plains  ecoregions.  In  the  plains, 
the  mean  overall  score  on  the  MTDEQ-RAM  was 
0.75  (B)  for  the  Northwestern  Great  Plains  and 
0.79  (B)  for  the  Northwestern  Glaciated  Plains.  In 
the  mountains,  the  mean  overall  score  was  0.92  (A) 
for  wetlands  in  the  Northern  Rockies,  0.93  (A)  for 
wetlands  in  the  Middle  Rockies,  and  0.93  (A)  for 
wetlands  in  the  Canadian  Rockies. 

We  evaluated  land  cover  and  land  use  metrics 
for  the  entire  data  set  to  identify  metrics  with  a 
significant  relationship  to  overall  wetland  condition 
score  and  an  R-squared  value  of  0.20  or  greater. 
Several  metrics  were  significant  at  one  or  more 
scales,  but  no  individual  metric  had  an  R-squared 
value  at  or  above  our  target  level.  To  determine 
whether  we  could  find  a  combinations  of  metrics 
with  some  predictive  value,  we  selected  metrics 
that  were  significant  across  at  least  three  of  the  four 
levels  tested,  and  evaluated  them  using  backward 
stepwise  regression. 

At  the  6th  code  HUC,  1  kilometer,  500  meter,  and 
200  meter  levels,  the  combined  metrics  of  percent 
forest  cover,  road  density,  and  number  of  stream 
road  crossings  had  the  strongest  predictive  value 
for  overall  score.  However,  the  R-squared  value 
for  the  three  metrics  combined  was  only  0.14  at  the 
6th  code  HUC  level.  When  measured  at  the  buffer 
level,  the  combination  of  metrics  had  an  even 
lower  R-squared  value.  At  1.000  meters,  the  R- 
squared  value  dropped  to  0.13.  at  500  meters  it  was 
0.11.  and  at  200  meters,  the  R-squared  value  was 
0.09.  and  road  density  was  no  longer  significantly 
correlated  with  overall  score. 


The  low  and  declining  R-squared  values  in  the 
initial  analysis,  coupled  with  the  ecoregional  skew 
in  the  dataset.  suggested  that  the  datasets  might 
be  better  analyzed  in  two  groups.  From  our  field 
work,  we  knew  that  road  density  and  road  stream 
crossings  are  generally  lower  in  the  mountain 
ecoregions.  while  percent  forest  is  generally  higher, 
and  so  it  appeared  likely  that  the  analysis  was 
detecting  ecoregional  differences  rather  than  the 
impact  of  land  uses.  Using  ArcGlS.  we  divided  the 
data  into  a  plains  set  (Northwestern  Great  Plains 
and  Northwestern  Glaciated  Plains)  and  a  mountain 
set  (Northern  Rockies.  Canadian  Rockies,  and 
Middle  Rockies),  and  reran  the  analysis. 

Initial  Analysis:  Step  2 

With  the  data  divided,  percent  forest  was  no 
longer  significant  at  any  scale.  For  wetlands  in 
the  mountain  ecoregions  (n=262),  road  density 
was  the  only  metric  that  was  significant  at  all 
levels,  although  the  R-squared  value  was  never 
higher  than  0.07.  At  the  1 .000  meter  buffer,  the 
percentage  of  crop  agriculture  was  also  significant 
(p  <  0.001  R-squared  =  0.08).  aUhough  it  had  no 
significance  at  other  buffer  distances.  In  the  plains 
ecoregions.  none  of  the  metrics  were  significant 
at  200  meters.  Percent  natural  grassland  and  road 
density  within  500  meters  were  both  significantly 
correlated  with  overall  score  but  had  very  low  R- 
squared  value  (0.02  and  0.01.  respectively).  At  the 
1.000  meter  buffer  scale,  only  the  number  of  stream 
road  crossings  was  significant.  No  metric  was 
significantly  correlated  to  overall  wetland  condition 
when  measured  at  the  6th  code  HUC  level  in  either 
the  mountain  dataset  or  the  plains  dataset. 

Because  no  human  disturbance  factor  appeared 
to  explain  the  unbalanced  distribution  of  scores 
among  ecoregions,  we  turned  to  environmental 
factors.  The  pattern  in  the  data  suggested  that 
some  combination  of  precipitation  and  elevation 
was  at  play,  so  we  reassessed  the  data  using  REAP 
as  a  variable.  REAP  was  highly  significant  when 
evaluated  against  the  entire  dataset.  (R-squared 
=  0.2 1 ).  When  we  tested  the  predictive  value 


ot  REAP  within  the  two  ecosystem  groupings, 
it  was  only  significant  in  the  mountain  regions 
(R-squared  =  0.13).  At  the  wetland  sites  in  the 
mountain  ecoregions,  REAP  ranged  from  43  to 
206  centimeters  (1 7.2  to  82.4  inches),  with  a  mean 
value  of  104  centimeters  (41 .6  inches).  By  contrast, 
values  at  wetland  sites  in  the  plains  ecoregions 
were  much  more  tightly  grouped,  with  a  range 
from  27  to  62  centimeters  ( 1 0.8  to  24.8  inches)  and 
a  mean  of  35  centimeters  (14  inches). 

CART  Analysis 

The  CART  analysis  was  run  with  a  proportional 
reduction  in  error  (PRE)  target  set  at  0.05.  We 
used  ail  predictor  variables,  evaluating  the  data 
with  metrics  measured  in  the  1 .000  meter  buffer. 
However,  for  the  dataset  as  a  whole,  no  factor  other 
than  REAP  reduced  the  PRE  an  additional  0.05. 
The  same  was  true  when  we  set  the  PRE  to  0.03. 
Dividing  the  data  into  mountain  and  plains  groups 
failed  to  improve  performance.  At  a  PRE  of  0.05, 
REAP  was  the  only  splitting  factor  in  mountain 
sites,  while  no  splitting  could  be  achieved  in 
the  plains  sites.  Because  we  were  searching  for 
predictor  variables  based  on  human  disturbance, 
we  then  removed  REAP  as  a  predictor  and  ran  the 
analysis  again.  In  the  mountains,  road  density 
equal  to  zero  emerged  as  the  first  splitting  factor, 
with  percent  natural  grassland  greater  than  2.5%  as 
the  second.  In  the  plains,  no  factor  emerged  until 
we  set  the  PRE  to  0.02,  where  number  of  stream 
crossings  produced  a  split. 

We  had  similar  results  with  all  other  buffer 
distances.  In  the  mountains,  the  splitting  factors 
remained  the  same.  In  the  plains,  stream  road 
crossings  dropped  out  entirely  at  500  meter  and 
200  meter  buffers.  No  factor  split  the  data  until 
road  density  emerged  when  we  dropped  the  PRE  to 
0.01. 


In  short,  the  CART  merely  reaffirmed  the 
observations  made  with  linear  analysis.  Most  of 
the  variation  in  our  dataset  was  not  explained  by 
measured  environmental  factors  or  human-driven 
landscape  factors,  suggesting  that  site-specific 
factors  may  have  the  most  pronounced  impact 
on  wetland  condition^.  However,  we  were  also 
cognizant  that  the  dataset  itself  was  not  ideal  for 
building  a  data-driven  model.  Because  of  its 
design,  many  of  the  datapoints  were  concentrated 
in  roadless  areas  with  no  measurable  human 
impacts.  Others  were  clustered  in  areas  where 
road  density  was  low  but  grazing  pressures  were 
high.  Despite  emerging  as  a  significant  predictor  of 
wetland  condition  in  the  mountain  ecoregions.  crop 
agriculture  occupied  10%  or  more  of  the  1,000 
meter  buffer  at  only  four  assessment  sites. 

The  Montana  Landscape  Integrity 
Model 

Our  Montana  Landscape  Integrity  Model  (MTLIM) 
is  an  inverse  weighted  distance  model  premised 
on  the  idea  that  ecosystem  processes  and  functions 
achieve  their  fullest  expression  in  areas  where 
human  activities  have  the  least  impact.  In  the  case 
of  wetlands,  it  presumes  that  wetland  condition  will 
be  poorest  in  close  proximity  to  roads,  commercial 
or  industrial  development,  urban  areas,  resource 
extraction  sites,  or  hydrologic  modifications. 

The  model  is  built  into  a  single  raster  layer 
covering  the  entire  state.  The  raster  pixel  size  is  30 
meters  by  30  meters,  or  900  square  meters.  Pixel 
values  range  from  100  to  745.  Pixels  with  a  value 
of  1 00  are  beyond  the  furthest  extent  of  any  buffer 
distance  for  any  impact  (Figure  2).  Conversely, 
pixels  with  high  values  (>500)  are  typically  in 
urban  areas  at  the  intersection  of  multiple  stressors 
(Figure  3). 


■*  We  made  several  other  attempts  to  extract  patterns  from  the  data  using  both  CART  and  linear  techniques.  For  example,  because 
the  scores  from  the  randomly-drawn  dataset  were  not  evenly  distributed  -almost  half  the  wetlands  were  ranked  "A"  with  a 
score  of  0.9  or  more — we  drew  a  geographically  stratified  sample  of  40  wetlands  in  each  rank  category  from  the  overall  dataset 
and  repeated  all  the  analyses  described  here.  Again,  only  the  factors  that  were  seen  as  significant  in  the  areas  and  at  the  scales 
described  above  stood  out.  and  R-squared  values  were  within  the  same  ranges.  We  also  added  additional  factors,  evaluating 
land  ownership,  septic  density,  population  density,  management  status,  and  water  rights  points  of  use.  Although  many  factors 
were  positively  correlated  to  wetland  condition,  no  factor,  alone  or  in  combination,  had  notable  predictive  value.  Those  having 
significance  at  several  scales  were  incorporated  into  the  expert  model. 


Figure  2.    Wetland  uiihiii  high- integrity  landscape. 


Figure  3.   iVetland  adjacent  to  multiple  stressors. 


We  calculated  Landscape  Integrity  Model  (LIM) 
scores  at  three  buffer  distances,  and  regressed 
these  scores  against  the  raw  MTDEQ-RAM  overall 
scores.  The  mean  LIM  score  calculated  from  pixels 
within  100  meters  of  the  assessment  point  yielded 
the  highest  R-squared  value,  0.08.  The  LIM  did 
not  account  for  significantly  more  variabilitv'  than 
any  of  the  individual  metrics  or  combinations  from 
the  original  analysis.  Again,  predictive  value  was 
ver\  poor  in  the  plains  ecoregions  (R-squared  = 
<0.01 ).  In  the  mountain  ecoregions,  the  predictive 
value  of  the  LIM  was  stronger  (R-squared  =  0.21). 

Using  CART  analysis,  we  identified  the 
combinations  of  REAP  and  LIM  scores  that  best 
accounted  for  the  variability  in  the  data.  After 
several  iterations,  we  found  that  the  formula 
yielding  the  best  results  was  as  follows  (all  REAP 
values  are  in  centimeters): 


A  =  (Mean  LIM  <120  and  REAP  >30)  OR  (Mean 
LIM  <180  and  REAP  >50)  OR  (Mean  LIM  <220 
and  REAP  >60) 

B  =  (Mean  LIM  <180  and  >120  and  REAP  <50) 
OR  (Mean  LIM  <i20  and  REAP  <30)  OR  (Mean 
LIM  <120  and  >180  and  REAP  <60  and  >40)  OR 
( Mean  LIM  <400  and  >220  and  REAP  >30) 

C  =  (Mean  LIM  <  220  and  >  i  80  and  REAP  <40) 
OR  (Mean  <350  and  >220  and  REAP  <30)  OR 
(Mean  LIM  >400  and  REAP  >50) 

D  =  Mean  LIM  >400  and  REAP  <50 

When  we  compared  the  rank  assigned  w  ith  these 
formulae  to  the  rank  assigned  with  the  MTDEQ  - 
RAM  on  the  original  591  data  points  and  test  datset 
of  1 00  points.  49.6%  were  classified  accurately. 
Accuracy  in  the  mountain  ecoregions  was  higher 
(74.8%).  When  applied  to  the  180  validation 
dataset  points,  5 1  %  were  classified  correctly 
across  all  ecoregions.  and  64.5%  were  classified 
accurately  in  the  mountains.  With  both  the 
calibration/test  and  validation  datasets.  accuracy 
was  highest  when  predicting  A-ranked  wetlands 
in  mountain  areas  (94%  with  the  calibration/test 
data  and  75%  with  the  validation  data).  In  the 
validation  dataset,  65%  of  the  plains  wetlands  and 
83.5%  of  the  mountain  wetlands  were  classified 
acceptably,  (i.e..  the  classifications  were  either 
accurate  or  no  more  than  one  rank  higher  than  what 
was  assigned  in  the  field) 

When  compared  to  the  Ecological  Integrity 
Assessment  scores,  the  classifications  were 
accurate  in  60  of  108  cases  (55.5%)).  even  though 
most  of  these  assessments  were  conducted  in  the 
plains  ecoregions.  74.07%  of  the  assigned  ranks 
were  acceptable.  The  best  performance  was  with 
B-ranked  wetlands,  where  45  of  70  (64.2%)  were 
correctly  classified  and  88.5%  were  classified 
acceptably. 

In  the  dry-run  Level  I  assessment  of  1,800 
wetlands  that  had  been  randomly  selected  for 
this  project.  45%  were  ranked  as  A,  44%  were 
ranked  as  B.  7%  were  ranked  as  C  and  4% 
were  D.  When  only  mountain  ecoregions  were 
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considered  (n=l,187),  56%  were  ranked  as  A,  35% 
were  ranked  as  B,  3%  were  ranked  as  C,  and  6% 
were  ranked  as  D.    This  regional  distribution  of 
predicted  ranks  differs  from  the  distribution  of 
measured  scores  for  the  591  points  used  in  the 
analysis,  mostly  in  the  B  and  C  ranks.  Of  the  591 
datapoints  assigned  ranks  based  on  MTDEQ-RAM 
scores,  48%)  were  A,  25%)  were  B,  21%)  were  C, 
and  6%  were  D.    In  the  mountain  ecoregions, 
75%  were  A,  16%  were  B,  8%  were  C,  and  less 
than  I  %  were  D.  However,  we  note  that  the 
stratification  schemes  underlying  the  two  datasets 
were  very  different.  The  amphibian  survey  dataset 


from  which  the  591  samples  were  drawn  explicitly 
stratified  along  ownership  gradients,  and  a  high 
percentage  of  the  datapoints  in  the  mountains  were 
drawn  from  managed  public  lands.  By  contrast,  the 
sampling  scheme  for  the  1,800  wetlands  classified 
in  this  study  did  not  use  ownership  as  a  sfrata.  It 
may  be  that  the  higher  percentage  of  A-ranked 
mountain  wetlands  in  the  MTDEQ-RAM  dataset 
is  explained  by  the  fact  that  most  occurred  on 
public  land.  The  wetlands  selected  for  this  project 
occurred  across  all  ownership  strata,  and  included 
many  wetlands  in  proximity  to  developed  mountain 
valleys  in  western  Montana. 


12 


Discussion 


This  study  demonstrated  the  potential  of  landscape- 
level  metrics  and  models  to  predict  wetland 
condition  using  remotely-sensed  data  in  Montana. 
At  the  same  time,  it  showed  that  environmental 
variables  and  site-specific  activities  may  be  far 
more  important  drivers  of  wetland  condition  than 
land  uses  occurrina  at  a  broader  scale. 


metrics  involving  distance  from  crop  agriculture 
or  roads,  the  two  land  uses  where  soil  disturbance 
is  most  common,  could  be  rescaled.  This  would 
be  most  successful  if  it  was  based  on  additional 
field  data  that  could  support  development  of  a  clear 
relationship  between  buffer  distance  and  presence 
of  exotics  and  invasive  plants  in  wetlands. 


Further  development  of  a  landscape- level 
assessment  tool  will  require  a  more  robust  set 
of  field  data.  The  MTDEQ-RAM,  while  useful 
for  quick  characterization  of  wetland  impacts, 
is  problematic  insofar  as  it  conflates  stressors 
and  condition.  For  example,  across  the  entire 
state,  proximity  to  roads  and  road  density  within 
buffers  were  the  only  landscape  metrics  that  had 
a  consistent  negative  correlation  with  overall 
wetland  condition  measured  with  the  MTDEQ- 
RAM.  However,  proximity  to  roads  is  a  metric  in 
the  MTDEQ-RAM  buffer  condition  score,  so  the 
correlation  may  rest  in  large  part  on  near-identical 
metrics  rather  than  on  any  real  causal  relationship 
between  road  density  and  wetland  condition.  As 
we  build  new  datasets  using  rapid  and  intensive 
Ecological  Integrity  Assessments,  which  separate 
condition  metrics  and  human  disturbance  metrics, 
it  will  be  useful  to  revisit  this  analysis.  The  fact 
that  the  Landscape  Integrity  Model  and  REAP  data 
were  better  able  to  predict  EIA  scores  in  the  plains 
ecoregions  than  MTDEQ-RAM  scores  is  promising 
in  this  regard. 


We  did  not  do  metric-by-metric  correlations 
between  landscape- level  stressor  variables  and 
wetland  condition  variables.    The  MTDEQ-RAM 
metrics  are  semi-quantitative,  but  they  rely  heavily 
on  subjective  assessment  by  the  user.  For  example, 
question  3.2,  under  "Hydrogeomorphology 
Condition"  asks  the  user  to  assess  "[d]egree  of 
wetland  habitat  negatively  altered  by  addition 
or  withdrawal  for  irrigation,  livestock  watering, 
drainage,  etc."  with  the  possibilities  being  "Non- 
occuring/slight  (10  points),'"  "Moderate  (4  points)" 
or  "Severe  (0)."  Similarly,  the  "Vegetation 
Condition"  metric  4. 1  asks  the  user  to  evaluate 
"How  much  emergent  vegetation  is  impacted  by 
trampling  or  other  human-caused  disturbance?" 
with  the  choices  being  "None  present/minimal 
(10  points),"  "Some  present  6-15%  (8  points)," 
"Common  occurrence  16-25%  (4  points)"  or  "Very 
apparent  >25%  (0  points)."     If  the  assessment 
data  were  fully  quantified,  we  think  that  metric-to- 
metric  comparisons  would  help  to  sharpen  the  data 
and  to  tease  out  differences  between  ecoregions 
(Hychka  et  al.  2007). 


The  results  of  this  study  also  suggest  that  further 
model  development  should  separate  out  the  plains 
ecoregions  and  the  mountain  ecoregions.    In  the 
wide-open  landscapes  that  characterize  much 
of  the  plains,  airborne  weed  seeds  and  dust  can 
travel  much  further  than  in  forested  regions,  so 
the  distance  at  which  individual  landscape-level 
stressors  operate  may  be  greater.  In  our  field 
studies,  we  have  seen  several  cases  where  the 
presence  of  a  thickly  forested  or  shrubby  buffer 
between  a  weed  source  and  a  wetland  appears  to 
protect  the  wetland  from  the  spread  of  invasive 
plants.  By  contrast,  in  some  areas  of  the  plains, 
agricultural  weeds  or  exotics  are  present  in  nearly 
all  wetlands.  This  suggests  that  in  the  plains. 


Of  course,  the  best  results  will  be  obtained  when 
Level  I  assessment  metrics  are  calibrated  against 
Level  III  data  collected  by  expert  users.  The 
MTDEQ-RAM  was  designed  for  use  by  field 
crews  with  only  enough  botanical  knowledge  to 
identify  a  small  suite  of  exotic  plant  species  (reed 
canarygrass,  smooth  brome,  quack  grass,  Kentucky 
bluegrass,  creeping  bent  grass,  meadow  foxtail, 
tall  fescue,  timothy,  sweet  clover,  and  Russian 
olive)  and  noxious  weeds  (tamarisk,  Canada  thistle, 
whitetop  cress,  spotted  knapweed,  leafy  spurge, 
purple  loosestrife,  yellowflag  iris,  and  Eurasian 
milfoil).    We  think  that  one  of  the  reasons  the 
Landscape  Integrity  Model  and  REAP  data  were 
more  accurate  when  predicting  ranks  obtained  with 
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EIAs  is  that  the  EIAs  require  users  to  be  able  to 
evaluate  the  degree  to  which  the  vegetation  departs 
from  reference  standard  communities  expected  in  a 
given  ecological  system,  and  to  be  able  to  quantify 
the  percent  of  all  exotics  and  invasive  plant  species 
present.  Although  EIAs  are  also  only  semi- 
quantitative, this  improvement  in  performance 
suggested  that  better  field  data  will  yield  better 
Level  I  assessment  metrics. 

We  should  note,  too,  that  we  calibrated  Landscape 
Integrity  Model  scores  to  yield  the  highest  overall 
accuracy,  but  that  our  results  may  favor  wetlands 
in  mountain  ecoregions  simply  because  that  is 
where  most  of  the  A-ranked  wetlands  in  the  dataset 
occur.  Again,  had  we  been  working  with  more 
precise  data,  it  would  have  been  worthwhile  to 
attempt  separate  calibrations.  However,  without 
statistically  significant  relationships  between 
wetland  condition  in  the  plains  and  most  of  the 
individual  landscape  metrics  that  went  into  the 
LIM,  this  seemed  premature. 

We  also  have  not  found  a  way  to  capture  the 
effects  of  grazing,  which  can  be  both  a  site-specific 
stressor,  as  when  livestock  enter  the  wetland,  or 
a  landscape- level  stressor,  when  they  remove 
vegetation  fi-om  the  buffer  area.  Currently,  no 
large-scale  GIS  grazing  layer  exists.  On  a  local 
scale,  it  might  be  possible  to  characterize  grazing 
intensity  using  animal  unit  months  (AUMs) 
obtained  from  land  management  agencies  if 
the  area  in  question  is  under  their  jurisdiction. 
Similarly,  it  may  be  possible  to  use  image 
classification  to  parse  out  rangelands  into  condition 
classes  based  on  spectral  signatures,  but  this  has 
not  been  done  for  the  state  as  a  whole,  and  is  a 
time-consuming  and  expensive  effort  that  itself 
requires  substantial  ground  truth  data. 

One  of  the  questions  raised  by  this  study  is  whether 
plains  wetlands,  in  general,  are  in  poorer  condition 
than  mountain  wetlands.  This  was  suggested  both 
by  the  MTDEQ  data  and  by  running  the  model  on 
a  randomly  selected,  spatially  distributed  set  of 


wetlands  across  the  state.  We  note  several  reasons 
that  this  may  be  so.  First,  in  a  period  of  protracted 
drought,  wetlands  in  more  arid  regions  would  be 
more  severely  impacted.  Wetland  obligate  plants 
will  be  more  stressed,  and  competition  with  both 
facultative  wetland  and  facultative  upland  species 
will  increase.  In  cases  where  wetland  obligate 
species  have  low  palatability  to  livestock  and 
wildlife  and  facultative  upland  plant  palatability  is 
high,  drought-induced  shifts  in  plant  communities 
(e.g.,  a  shift  from  an  Eleocharis  palustris 
association  to  a  Pascopyrum  smithii-Eleocharis 
spp.  association)  may  increase  grazing  pressure.  A 
second  explanation  lies  in  the  simple  fact  that  wet 
places  in  dry  areas  attract  more  human  disturbance 
in  the  form  of  water  impoundments,  diversions 
and  stock  water  development.  Although  we 
selected  only  wetlands  that  field  crews  had  labeled 
"natural"  for  our  analysis,  examination  of  aerial 
photography  and  NWI  maps  indicated  that  many 
of  these  "natural"  wetlands  were  dammed,  diked, 
or  receiving  water  flow  from  upstream  reservoirs. 
A  third  factor,  related  to  the  first  two,  is  resilience. 
Drought-stressed  wetlands,  or  seasonally  flooded 
wetlands  undergoing  earlier  drawdown  will  have 
more  bare  ground  for  longer  periods,  facilitating 
establishment  of  exotics  and  invasive  plant 
species.  The  same  is  true  of  wetlands  disturbed 
by  direct  human  influences.    Even  if  livestock 
is  removed  or  hydrological  modifications  are 
abandoned,  it  may  be  some  time  before  soil  and 
moisture  conditions  allow  for  a  return  of  native 
wetland  plant  communities.    Therefore,  we  think 
it  is  likely  that  plains  wetlands  are,  overall,  more 
likely  to  be  in  good-to-fair  condition  rather  than 
the  excellent-to-good  condition  found  in  the 
mountains.  However,  we  also  note  that  most  of  the 
intensive  human  developments  (cities,  recreation 
areas)  found  in  Montana  occur  in  the  mountains, 
and  that  a  different  sampling  scheme,  geared  to 
characterizing  wetland  condition  across  a  basin 
or  ecoregion  rather  than  to  surveying  amphibians, 
would  be  necessary  to  generate  a  dataset  that  could 
answer  this  question  conclusively. 
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How  TO  Use  the  Landscape  Integrity  Model 


Landscape-level  assessments  are  a  rough  surrogate 
for  site-specific  measurements  of  human-induced 
disturbances  that  may  atTect  wetland  condition 
(Hychka  et  al.  2007.  Wardrop  et  al.  2007).    Fhey 
are  not  themselves  assessments  of  condition  and 
cannot  substitute  for  site-specific  evaluations. 
However,  there  are  several  cases  where  use  of  the 
Landscape  Integrity  Model  and  REAP  data  would 
be  appropriate: 

1 .  Characterize  and  compare  wetland  disturbance 
or  likely  wetland  condition  across  a  geographic 
area.    Ranks  derived  from  the  LIM  and  REAP  data 
can  be  displayed  at  multiple  scales.  At  the  full  state 
level,  such  displays  can  illustrate  the  distribution 


of  stressors  across  regions  of  the  state.  Figure  4, 
for  example,  shows  the  1,800  randomly  selected 
points  scored  for  this  study.  A-ranked  wetlands  are 
clustered  in  the  mountain  regions,  while  B-  and 
C-ranked  wetlands  are  more  concentrated  in  the 
plains.  Similar  characterizations  could  also  be 
made  at  smaller  scales  if  digital  wetland  maps  are 
available'. 

2.  Characterize  watersheds  or  other  geographic 
units.  If  wetland  maps  are  not  available,  or  if  users 
simply  want  to  characterize  a  landscape  rather  than 
individual  wetlands,  then  the  LIM  can  be  used  as  a 
general  landscape  assessment  tool.  For  example, 
watershed  assessments  typically  include  some  sort 
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•  D-RANK 


Figure  4.  Stodel-ranked  wetlands. 


-   If  digital  wetland  maps  are  available,  we  would  recommend  using  a  revised  LIM  incorporating  an  "altered  wetland"  raster. 
We  did  not  do  this  because  the  goal  was  a  model  that  would  be  usable  statewide;  if  we  had  been  making  a  model  for  a  specific 
subbasin.  with  wetland  mapping,  we  would  have  substituted  altered  wetlands  for  other  hydrologic  alterations. 
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of  synthesis  of  stressors,  threats,  or  disturbances 
(e.g.,  Vance  and  Stagliano  2007,  2008).  When 
these  are  carried  out  at  the  basin,  subbasin,  or 
watershed  level,  nested  hydrologic  units  can  be 
compared  and  contrasted.  Figure  5,  which  shows 
the  Ruby  Valley  subbasin  in  Southwest  Montana, 
displays  average  LIM  scores  calculated  for  each 
subwatershed,  and  provides  a  quick  snapshot  of 


the  distribution  of  stressors  across  the  subbasin. 
More  precise  versions  of  the  display  can  show 
LIM  scores  for  each  pixel,  allowing  users  to 
focus  on  more  specific  areas.  Figure  6  depicts  the 
same  subbasins  with  non-aggregated  LIM  scores, 
allowing  users  to  visualize  areas  of  highest  and 
lowest  integrity. 
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Figure  5.  Landscape  integrity  by  subwatershed. 
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Figure  6.  LI M  values  by  pixel. 

3.  Identify  high-integrity  hot  spots  where  high 
quality  wetlands  are  likely  to  be  found.  Although 
site-specific  disturbances  will  always  be  a 
possibility,  the  LIM  can  be  used  to  locate  areas 
where  large-scale,  negative  human  influences  are 
minimal.  To  create  Figure  7,  we  used  a  conditional 
statement  in  Spatial  Analyst  to  identify  pixels 
with  a  LIM  value  of  100  (no  human  impacts). 
Then  we  used  block  statistics  in  Spatial  Analyst 


to  calculate  average  scores  over  3  kilometer  by 
3  kilometer  blocks.  Block  statistics  could  be 
calculated  for  larger  or  smaller  areas,  but  this  area 
allows  a  user  to  pick  out  contiguous  areas  of  high 
integrity.  Figure  8  is  zoomed  into  northwestern 
Montana,  where  the  populated  Flathead  Valley 
and  surrounded  timberlands  contrast  sharply  with 
the  pristine  sections  of  Glacier  National  Park  and 
the  Middle  Fork  Flathead  watershed.  In  this  case. 


17 


the  map  would  direct  searches  for  high-quality 
wetlands  (or  other  natural  features)  towards  the 
blocks  defined  as  highest  integrity. 

4.  Combine  with  other  tools  to  conduct  desk-based 
assessments  of  individual  wetlands.  The  wetland 
in  Figure  9  was  classified  to  an  A-rank  with  the 
LIM  score  and  REAP  data.  In  the  field,  it  was 
ranked  as  a  C.  Examination  of  the  high-resolution 
imagery  taken  in  2005  explains  why.  First,  the 
wetland  itself  has  been  excavated,  as  have  several 
in  the  area.  Cattle  paths  are  visible  to  the  west  of 


the  wetland,  and  a  road  (not  depicted  on  any  GIS 
layer)  is  evident  to  the  north.  If  LIM  and  REAP 
scores  are  used  as  individual  metrics  and  additional 
metrics  are  developed  (e.g.,  evidence  of  grazing 
within  a  1-km  buffer,  percent  altered  wetlands 
within  a  1  -km  buffer,  length  undigitized  roads  or 
two-tracks,  etc),  then  accuracy  would  be  greatly 
improved  in  many  cases.  The  same  would  be  true 
in  forested  areas,  where  clearcuts  and  fires  may 
have  occurred  after  the  GIS  datasets  were  created, 
but  prior  to  acquisition  of  the  high-resolution 
imagery  in  2005. 
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Figure  7.  High-integrity  areas. 
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Figure  8 .  Northwestern  Montana  LIM  blocks. 


Figure  9.  Misclassified  wetland. 
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Conclusions 


Human  disturbances  can  affect  wetlands  at  multiple 
scales,  and  because  they  tend  to  interact  both  with 
each  other  and  with  environmental  variables,  it  is 
exceedingly  difficult  to  quantify  or  predict  their 
effects.  In  many  parts  of  Montana,  environmental 
variables,  particularly  precipitation,  may  be  the 
most  influential  drivers  of  wetland  condition. 
Nevertheless,  by  integrating  human  disturbances 
through  a  GIS,  we  were  able  to  build  a  model  of 
landscape  integrity  that  can  be  combined  with 
precipitation  data  to  generally  categorize  wetlands 
into  condition  ranks.  While  we  would  caution 
against  using  this  approach  to  rank  any  one 
wetland,  we  think  that  it  offers  a  way  to  quickly 
characterize  the  condition  of  wetlands  across  an 
entire  subwatershed  or  watershed  and  to  compare 
wetland  condition  in  different  geographic  areas. 


Ultimately,  landscape-level  wetland  assessment 
may  never  be  a  satisfactory  substitute  for  field 
assessments.  Even  in  rural  states  like  Montana, 
pockets  of  rapid  development  soon  make  derived 
image  classification  products  like  the  NLCD 
obsolete.  Similarly,  natural  stressors  like  insect 
outbreaks  and  disease  are  having  widespread 
impacts  on  our  forests,  and  salvage  logging 
operations  are  likely  to  increase  in  coming  years. 
Nonetheless,  when  field  assessments  are  not 
practical,  landscape-level  assessment  is  the  only 
alternative.  With  completion  of  the  Landscape 
Integrity  Model  and  an  approach  for  integrating  it 
with  precipitation  data,  we  now  have  a  tool  offering 
a  quick  and  efficient  way  to  begin  that  task. 
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Appendix  A.  MTDEQ  Rapid  Assessment  Form 


Montana  DEQ  -  Wetland  Rapid  Assessment  Form  (version 20) 


Site  Number 


Assessment  Number 


Site  Name 


Date 


Land  Ownership 


Pefson(s)  Assessing  Wetland  &  Affiliations 


HUC  4»'/5»'  Code 


HUC  WS*-  Name 


Elevation  (ft) 


Location  Infonnation 


UTME 


UTMN 


Datum 


GPS  ID 


NAD27 
NAD83 
Other 


UTMZone 


11 
12 
13 


GPS  error  (include  units) 


General  Site  Description  (location .  WMWe  observations,  Beaver  Acti\rity,  Outstanding  Features,  Vegetative  Types,  observed  impacts,  etc.): 


Photos 

Photo  # 

Direction  Facing              Description  of  what  is  in  the  photo 

1.0  Wetland  Classification 


1.1  Wetland  Is  being  assessed  to  reflect  (Circle) 

1 .2  HGM  Classification  (Circle  one  Class  or  Subclass) 

Nalural  Wetland  Type  (assess  potential) 

Altered  Wetland  Type  (assess  capability) 

Completely  Altered  (no  longer  lunclioning  as  a  wetland, 
and  it  is  not  leasible  to  survey  wetland  condition) 

'What alterations  have  tieen  made? 

Rivenne 

Upper  Perennial 

Lower  Perennial 

Non-Perennial, 
Intermittent  or 
Ephemeral 

Depressional 

Closed 

Open  groundwater 

Open  surface  water 

Lacustrine  Fringe 

Slope 

Open  Spring 

Riverine  Spring 

Fen 

Wet  Meadow 

Mineral  Soil  Flats 

Playa  Lalies 

Appendix  A  -  1 


1.3  Cowardin  Wetland  Classification  {Note:  wetlands  sites  can  have  more  than  one  system) 

Identify  a  System,  Subsystem,  Class.  Water  Reqime,  Modifier  (if  present),  and  the  percent  cover  of  all  categories  present 

System 

Subsystem 

Class 

Water 
Regimes 

Modifiers 

Percent 

Oetermine  the  wetland  area 

by  locating  the  boundary 

where  wetland  dependent 

vegetation  meets  vegetation 

and  features  not 

characteristic  of  wetlands 

(See  guidebook  for  more 

information) 

Do  not  include  limnetic 
subsystems  which  are  deep 

water  habitats  that  are 

greater  than  2  meters  (6.6 

feet)  or  the  maximum  extent 

of  nonpersistent  emergents. 

It  these  grow  at  depths 

greater  than  2  m. 

Types  of  Water  Regimes  and  Modifiers 

Riverine 
(SUcamj 

Lower  Perennial 
(Larger  Tnbutary) 

Rodcy  Bottom 

Water  Regimes  -  Choose  the  regime  that 

is  most  common  in  the  area. 

A  Temporarily  Flooded 

B  Saturated 

C  Seasonally  Flooded 

D  Seasonally  Flooded/Well  Drained 

E  Seasonally  Flooded/Saturated 

F  Semipermanently  Flooded 

U  Unknown 

Modifiers 

g  excavated 
h  impounded 
i  diked 

1  partly  drained 
k  fanned 
1  artificial  dam 
m  beaver  dam 
0  diverted 
p  rip  rap 

Unconsolidated  Bottom 

Aquatic  Bed 

Emeijent  Wetland 

Rocky  Shore 

Unconsolidated  Siiore 

Upper  Perennial 
(Smaller  Tributary) 

Rodty  Bottom 

Unconsolidated  Bottom 

Aquatic  Bed 

RodiyStioTB 

Unconsolidated  Shore 

Intermittent 

Stream  Bed 

Lacustrine 
(Lake) 

Limnetic 
(Deepwalei  habilaO 

Rocky  Bottom 

Unconsolidated  Bottom 

Aquatic  Bed 

Littoral 

(Between  Shore  and 

Dccpwater  Habrlal) 

Rocky  Bottom 

Unconsolidated  Bottom 

Aquatic  Bed 

Rod<y  Shore 

Unconsolidated  Shore 

Palustrine 
(Pond  or  riparian) 

Rocky  Bottom 

Aquatic  Bed  =  plants  growing  in  water 

Rocky  Bottom/  Shore  =  cobble  or  rock 

along  Shore 

Unconsolidated  Bottom/  Shore  =  muddy 

Emergent  =  grasses,  sedges,  nishes.  etc. 

Scrijb-Shn^b  =  Bushes,  Vegetaton  less 

than  20ft  tall 

Forested  =  woody  vegetation  that  is  6  m 

tall  or  taller 

Unconsolidated  Bottom 

Aquatic  Bed 

Emetjent  Wetland 

Rocky  Shore 

Unconsolidated  Shore 

Moss-Lichen  Wetland 

ScnjtvShnjb  Wetland 

Forested  Wetland 

2.0  Site  Characterization 

2.1  Are  Fish  Present?    |  Yes  |         1  No    |         |NotSure       |         [Species  (if  known)?! 

2.2  Amphibian  and  Aquatic  Reptile  Species  Observed  -  check  and  describe  life  stage  below:  Eqqs,  tadpole,  adult 

Common  Name 

Life  Stage 

Common  Name 

Life  Stage 

Common  Name 

Life  Stage 

Boreal  Chonis  Frog 

Snapping  Turtle 

Long-toed  Salamander 

Bullfioq 

Spiny  Softshell 

Northem  Leopard  Frog 

Coeur  D'Alene  Salamander 

Tiger  Salamander 

Pacific  Treefrog 

Columbia  Spotted  Frog 

Western  Hognose  Snake 

Painted  Turtle 

Common  Gartersnake 

Tenestrial  Gartersnake 

Plains  Garter  Snake 

Great  Plains  Toad 

Western  Toad 

Plains  Spadefoot 

Western  Skink 

Woodhouse's  Toad 

Rocky  Mtn  Tailed  Frog 

Smooth  Greensnake 

Other  (describe  if  unknown):                                                                                                                                   1 

2.3  Estimate  the  Percent  of  Standing  Water                                                                                                                                               | 

Percentage  of  standing  water  body  <  50  cm  depth 

0 

1-25 

26-50 

51-75 

76-100 

Percentage  of  standing  water  body  50-200  cm  depth 

0 

1-25 

26-50 

51-75 

76-100 

Percentage  of  standing  water  body  >200  cm  depth 

0 

1-25 

26-50 

51-75 

76-100 

2.4  Threatened  or  Endanc 

lered  Species  Observed  -  check  if  present  and  describ 

e  in  the  space  provided  below 

Check 

Species 

Region  Found 

Status 

Least  Tem 

Near  Fori  Peck  Dam  &  Miles  City 

Endangered 

Whooping  Crane 

Northeastern  Montana 

Endangered 

Bald  Eagle 

Entire  region 

Threatened 

Piping  Plover 

North-central  and  Eastern  portions  of  the  state 

Threatened 

Black-Footed  Ferret 

Northeastern  Montana 

Endangered 

Canada  Lynx 

Entire  region 

Threatened 

Gray  Wolf 

Entire  region 

Threatened/Endangered 

Grizzly  Bear 

Greater  Yellowstone,  Northern  Continental  Divide.  Cabinet- Yaak.  Bitte 

iToot  Selway  Ecosystems 

Threatened 

Bull  Trout 

Entire  Region 

Threatened 

Pallid  Sturtieon 

Fort  Peck  S  Yellowstone  River  below  Powder  River  mouth 

Endangered 

White  Sturgeon 

Kootenai  River 

Endangered 

Water  Howellia 

Northwestern  Montana 

Threatened 

Ule  Ladies'  -Tresses 

Southwest  and  Southcentral  Montana 

Threatened 

Please  comment  on  what  was  observed  (scat,  tracks,  etc.):            | 
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2.5  Check  ami  of  surface  area  of  any  emergent  vegetation 


Type 


Sedges 


Cattails 


Grasses 


Rushes 


Waterlilies 


Shrubs 


Trees 


1-25% 


25-50% 


50-75% 


76-100% 


'////////f///////f4 


Q 


Grasses 

Sedges 

Rushes 
Fence 


-0 


O 


Irees 
I'lioto 

Shrubs 

Assessment 
Boundaiy 


Other 


Please  describe; 


2.6  Site  Map  for  Wetland  Assessment  Area 

(site  map  can  be  substituted  with  a  high-resolution  aerial  photo) 


For  Riverine  sites:  include  !ength=  1 00m,  width=as  wide  as  outermost  meander.  For  all  other  sites;  1 00  m  ■  1 00m  or  the  entire  wetland,  if  smaller. 
Buffer  occupies  100m  on  either  side  of  the  wetland.  Specifics  for  determining  assessment  area  are  available  in  the  handbook. 
Grid  Scale:  1  square  =           m 

N 

■J 

•      Note  all  photo  locations  and  directions  What  is  the  overall  size  of  the  wetland? 

X 
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3.0  Hydrogeomorphology  Condition 


Degree  of  hydrologic  disturbance 
(All  Wetland  Types) 


Non  Occurring/Slight 


Moderate 


Severe 


3.1  Degree  of  wetland  surface  or  subsurface  flow  patterns  that  has 
been  "negatively'  altered  by  human  disturtance  (e.g..  roads, 
buildings,  rip  rap,  levees,  bridges  approaches,  weirs,  dams,  etc) 

•Consider  how  structures  accommodate  safe  passage  of  flows 
(e.g„  lower  the  ratng  it  headcuts  are  affecting  dam  or  spillway) 


10 


3.2  Degree  of  wetland  habitat  negatively  altered  by  addition  or 
withdrawal  for  irrigation.  livestoci<  watering,  drainage,  etc 

'Consider  impacts  from  any  abnormal  fluctuating  water  levels 


10 


3,3  Amount  of  wefland  habitat  negatively  alleted  by  dredging  or 
fliling 


10 


3,4  Percent  of  assessment  area  and  the  degree  to  which  the 
wetland  is  disturbed  by  pugging  or  hummoci<ing  from  animal 
hooves 

Slight=  Pugging  is  minimal  or  shallowWummocidng  has 
occun'edA/egetation  and  banl(  stability  is  intact  or  recovering 
Moderate^  Pugging  is  minimal/Hummoclis  are  deep/Wetland  is 
beginning  to  dry  out  Severe=  Hummocks  are  deep/  Pugging  is 
commonA/egetation  is  dead  or  absent 


<=25% 

None  Occurring  10 

Sliglit  9 

Moderate  6 

Severe  5 


26-75% 


76-100% 


Slight 

7 

Slight 

5 

Moderate 

4 

Moderate 

3 

Severe 

2 

Severe 

1 

Hydrogeomorphic  Condition  Index 

For  hydrologic  disturbance  take  the  sum  of  the  lowest  2  scores  (3,1-3.4)  and  divide  by  20: 


'Riverine  Index 


/20=( 


*For  Riverine  Sites  use  average  of  Riverine  and  Hydrogeomorphology  Indexes. 

Please  provide  comments  for  any  impacts  that  scores  <  5: 


Appendix  A-  4 


Hydrogeomorphology  -  Riverine  Wetland  Addendum  (Include  only  for  Riverine  Wetlands) 


The  at/t/j/ score  reflects  current  condition,  and  the  pvlMlidl\s  the  score  that  reflects  the  site  without  human  disturbance  (usually  the  maximum  score 

)■ 

3.5  Riverine  -Downcutting/lncisement:  Note:  The  presence  of  active  headcuts  should  nearly  always  keep  the 
stream  reach  from  t)einq  rated  sustainable. 

Actual 

Potential 

Stable  Channel 

6 

8 

Evidence  of  downcutting  ttial  is  beginning  to  stabilize 

6 

6 

Small  headcuts,  channel  is  in  beginning  staged  o(  unfavellng. 

4 

4 

Unstable  channel  that  is  incised  and  actively  widening,  banks  failure  is  common 

2 

2 

Deeply  tncised  resembling  a  gully 

0 

0 

3.6  Riverine  •  Percent  of  Stream  banks  iwith  active  Lateral  cutting: 

Actual 

Potential 

Lateral  bank  eroskjn  Is  in  balance  wilti  the  stream  and  its  setting 

8 

8 

There  is  a  minimal  amount  of  human-induced,  active  lateral  bank  erosion  occumng,  primarily  limited  to  outside  banks. 

5 

5 

There  Is  a  moderate  amount  of  human-induced  active  lateral  bank  erosion  on  either  or  both  outside  or  inside  banks 

3 

3 

There  is  extensive  humannnduced  lateral  bank  erosion  occumng  on  outside  and  inside  banks  and  straight  sections. 

0 

0 

3.7  Riverine  -  Stream  In  Balance  with  Water  and  Sediment  Supply:  Note:  Rosgen  B  and  naturally  occurring  D 
channels  are  exceptions. 

Actual 

Potential 

No  evidence  ol  excessive  sediment  removal  or  deposition,  or  that  the  stream  is  getting  wider. 

6 

6 

The  stream  has  widened  and/or  become  shallower  due  to  unstable  banks  or  from  de-watenng.  New  point  bars  are  often  forming  witfi  sitt  and 
sand  common 

4 

4 

The  stream  tends  to  be  very  wide  and  shallow.  Mid  channel  bars  are  often  presenl  (See  guidebook  for  prairies  streams  characteristics) 

2 

2 

The  stream  has  poor  sediment  transport  The  channel  is  often  braided  with  at  least  3  active  channels 

0 

0 

3.8  Riverine  -  Floodplain  Characterization:  (Rosgen  diagrams  are  available  in  the  handbook) 

Actual 

Potential 

Little  evidence  of  floodplain  erosion 

8 

8 

Floodplain  erosion  not  extensive 

6 

6 

Considerable  evidence  of  floodplain  erosion  and  occasional  headcuts 

4 

4 

Erosion  and  headcuts  within  the  floodplain  are  extensive.  Some  human<aused  stream  bank  erosion  is  occumng 

2 

2 

The  floodplain  is  very  limited  or  does  not  exist 

0 

0 

3.9  Riverine  •  Streambank  with  Vegetation  (Kind)  having  a  Deep,  Binding  Rootmass:  (see  Appendix  for 
stability  ratings  for  most  riparian,  and  other,  species) 

Actual 

Potential 

The  streambank  vegetaove  communities  are  compnsed  of  at  least  lour  plant  species  with  deep  binding  root  masses 

6 

6 

The  streambank  vegetative  communities  are  compnsed  of  at  least  three  plant  species  with  deep  binding  root  masses 

4 

4 

The  streambank  vegetative  communities  are  compnsed  of  at  least  two  plant  species  with  deep  binding  root  masses 

2 

2 

The  streambank  vegetative  communities  are  compnsed  of  one  or  no  plant  species  with  deep  binding  root  masses 

0 

0 

3.10  Riverine  -  Streambank  with  Vegetation  (Amount)  having  a  Deep,  Binding  Rootmass:  (see  Appendix  for 
stability  ratings  for  most  riparian,  and  other,  species) 

Actual 

Potential 

More  than  85%  of  the  floodplain  has  vegetation  with  a  stability  raCng  greater  than  or  equal  to  6 

6 

6 

75-  85%  of  the  floodplain  has  vegetation  with  a  stability  rating  greater  than  or  equal  to  6 

4 

4 

65-75%  of  the  floodplain  has  vegetation  with  a  stability  rating  greater  than  or  equal  to  6 

2 

2 

<  65%  of  the  floodplain  has  vegetation  with  a  stabil^y  rating  greater  than  or  equal  to  6 

0 

0 

Please  provide  comment  for  any  indiviiiual  score  <6: 

If  the  potential  is  not  at  maximum,  please  explain: 

Riverine  Index: 

1 

Sum  the  actual  scores  (3.5-3.1 0)  and  divide  by  the  sum  of  the  potential  scores  (usually  the  maximum  scores): 

Actual;              +         +           +           +          = 

* 

- 

Potential:          +           +         -i-           +          = 
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4.0  Veqetation  Condition  "Vegetation  should  only  be  assessed  within  the  wetland  assessment  area 

4.1  Bare  Ground 

None  present/  Minimal 
<=5% 

Some  Present 
6-15% 

Common  Occurrence 
16-25% 

Very  apparent 
>25% 

How  much  emergeni  vegetation  is 
impacted  by  trampling  or  other 
human-caused  disturtance? 

10 

8 

4 

0 

*For  Noxious  antj  Disturbance  Caused  Undesirable  plants,  look  to  the  abuntJance  of  harmful  species. 

4.2  Invasive  and  Disturbance  caused 

Some  small  patches  are 

Patches  are  large  or 

Patches  are  large  and 

undesirable  plants 

None  present 

often  present 

commonly  present 

extensive  or  Wetland  is 

(Rank  3  most  common  and  check  all  other 

<=5% 

6-25% 

Dominated 

observations) 

>25% 

Reed  Cana/y  grass                      Meadow  Foxtail 

Smooth  brome                          Tal  Fescue 

Quack  grass                             Timothy 

10 

7 

5 

2 

Kentucky  bkjeqrass                      Sweet  Ctover 

Creeping  Bent  grass                    Russian  Olive 

4.3  Noxious  Weeds 

Some  small  patches  are 

Patches  are  large  or 

Patches  are  large  and 

(Rank  3  most  common  and  check  all  other 

None  present 

often  present 

commonly  present 

extensive  or  Wetland  is 

observations) 

<=5% 

6-25% 

Dominated 

>25% 

Tamarisk  (Salt  Cedar) 

Leafy  Spuige 

Canada  Thistle 

Purple  Loosestrife 

10 

6 

3 

0 

White  Top  Cress 

YeHowflag  Ins 

Spotted  Knapweed 

Eurasian  Mtfoil 

Is  woody  vegetation  present?  Yes No "Sklp  the  rest  of  this  section  if  the  site  does  not  have  the  potential  for  tall  shrubs 

or  trees  or  woody  vegetation  isnotpresentdueto  natural  causes  (not  human  impacts  or  removal). 


4.4  Woody  Species  Establishment  and  Regeneration 


Actual 


Potential 


All  age  classes  of  desirable  woody  species  present  (see  Guidebook). 


10 


10 


One  age  class  of  desirable  woody  species  is  clearly  absent,  all  others  well  represented.  Often,  it  will  be  the  middle  age  group(s)  absent 


Two  age  classes  (seedlings  and  saplings)  of  native  shrubs  and/or  two  age  classes  of  native  trees  are  cleariy  absent,  or  ttie  stand  is  compnsed  of 
mainly  mature  species.  Other  age  classes  well  represented. 


DisturUance  induced,  (i.e.,  facultative,  facultative  upland  species  such  as  rose,  or  snowberry)  or  non-wetlands  dominate.  Woody  species  present 
consist  of  decadenVdying  individuals 


A  few  woody  species  are  present  (<  10%  canopy  cover),  but  herbaceous  species  dominate  (at  this  point,  the  site  potential  should  be  re-evaluated  to 
ensure  that  it  has  potental  tor  woody  vegetation).  OR,  the  site  has  at  2  5%  canopy  cover  of  Russian  olive  and/or  salt  cedar. 


0 


0 


4.5  Utilization  of  trees  and  shrubs: 


Actual 


Potential 


Few  to  none  of  the  available  second  year  and  older  stems  are  browsed 


10 


10 


Second  year  and  older  stems  lightly  browsed 


Second  year  and  older  stems  are  moderately  browsed. 


Second  year  and  older  stems  are  heavily  browsed.  Many  of  the  shrubs  have  either  a  "clubbed"  growth  fomi,  or  they  are  high-lined  or  umbrella  shaped. 


There  is  noticeable  use  (10%  or  more)  of  unpalatable  and  normally  unused  woody  species 


0 


4.6  Percent  of  physical  removal  of  tree/shrub  layer  or 
dead  wood  caused  by  concentrated  livestock  trampling  and  rubbing, 
drying  out  of  site  due  to  stream  incisemenL  human-caused  wetland 
drainage  or  flooding,  etc. 


<=5% 


10 


6-25% 


26-50% 


51-75% 


76-100% 


Please  provide  comments  for  any  individual  scores  less  than  6: 


If  Potential  is  not  at  maximum,  please  explain: 


Vegetation  Condition  Index 

Sum  all  scores  and  divide  by  the  total  possible  for  the  assessment  area.  60  for  sites  with  t/voody  species  (shrubs  and  tree);  30  for  sites  with  only  herbaceous 

vegetation). 

Only  Herbaceous  (4.1-4.3): + + = /30 


For  Herbaceous  and  woody  vegetation  (4.1-  4.6): 

(  /10+         /10+         /10  +  actual/potentialt  actual/potentlaU  /10)/6    = 
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5.0  Water  Quality:  Is  water  present?  Yes No *Skip  this  section  if  water  is  not 


present 


5.1  Alqae  and  Duckweed 

Latge  patches  means  50% 


Algae  growth  IS 
minimal 

10 


Algae  growtti  in  smail 
patches 


Algae  growth  in  large  patches 
4 


High  level  ol  algae  growth  in  continuous 
mats  with  odor  from  rotting  vegetation 

0 


S.2  Is  Wetland  Dominated  by  Cattails? 

•Dominated  means  70S 

Do  not  include  any  open  water  componenl 


Yes 


No       10 


5.3  Sediment  and  Turbidity 


5.3a  Is  there  evidence  of  excessive 
sediment  levels  caused  by  human 
activities?  (e.g.  bare  ground,  row  crops, 
erosion,  etc.  Do  not  include  trapped  sediment 
due  to  lieaver  damming) 


5.3b  Is  the  Water  Turbid? 


No  evidence  / 
Slight 

10 


No  Turtidity/ 
Slight 

W 


Moderate 
4 


Moderate 
8 


High 


High 


Average  Sediment  and  Turbidity  Score: 

- K= 


10     987654320 


5.4  Surface  oils  &  foams 

"Do  not  consider  sheen  (or  vegetation  decomposition 
(Should  be  evidence  ol  human  caused  source) 


No  evidence  ot  surface  oils 
or  foams 

10 


Evidence  of  surface  oils  or  foams 
3 


The  wetland  is  covered  with  surface  oils  or  (oams 
0 


5.5  Toxics-  (e.g.  Metals  from  mine  taiSngs, 
hydrocartjon  organic  materials,  or  Pesticides) 


No  evidence  of  tonics 


10 


Evidence  of  toxics,  however  aquatic  fife  is 
abundant  and  diverse 
5 


Evidence  of  toxics. 
Only  tolerant  aquatic  life  are  found 
0 


5.6  Salinity 

•Conductivity  measurements  are  not  necessary 


5.7  Are  saline  seeps  fallow  croplands,  oil  t>rines.  or 
severe  overgrazing  present  within  3  mites' 
Yes     No     Not  Sure 


No  evidence  of  saline  seeps 
Conductivity 
<  3000  uS/cm 
10 


Moderate  evidence  o(  saline  seeps 

Conductivity 

3000-15000  uS/cm 

5 


Significant  evidence  of  saline  seeps 

Conductivity 

>15000uS/cm 

0 


Water  Quality  Condition  Index:  Sum  the  lowest  2  scores  (5.1-5.6)  and  divide  by  20: 


/20: 


Please  comment  on  any  individual  scores  <  6: 


6.0  Buffer  Condition/  Deo 

ree  of  Stress 

stressors  in  100  meter  buffer 

None  present 
Very  few  present 
/Minimal 
Small  Patches 

Common 

Occurrence 

Large  patches  within 

Buffer 

Very  apparent  and  extensive 

Distribution 

Extensive  Large  Patches  throughout  entire  Buffer 

6.1  Amount  of  bare  ground 

10 

Slope 

Flat           6 
Moderate    4 
Steep          3 

Slope 

Flat            4 
Moderate    2 
Sleep          1 

Slope 

Flat=  <2  percent  grade 

Moderate=  2-10  percent  Grade 
Steep=  >10  percent  gradg 

6.2  Noxious  weeds 

(Use  Montana  Noxious  Weed  Pamphlet) 

10 

2 

0 

6.3  Disturbance-  caused  undesirable 
plants 

10 

4 

0 

Degree  of  Stress  in  Buffer 

None 

Occurrinq/Sliqht 

Moderate 

Severe 

6.4  Grazing  intensity 
in  100  meter  buffer 

10 

Slope 

Flat           7 
Moderate    5 
Steep          4 

Slope 

Flat           4 
Moderate    2 
Steep          1 

6.5  Recreational  Activities  (e.g. 
campground,  fishing  access  point, 
etc.) 

10 

Slope 

Flat            7 
Moderate    5 
Steep         4 

Slope 

Flat            4 
Moderate    2 
Steep         1 
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Percent  of  100m  buffer  occupied 
by  stressor 

0% 

1-25% 

26-50% 

>50% 

6.6  Hayfield 

10 

8 

6 

4 

6.7  Row  Crops 

10 

Slope 

Flat         7 
Moderate  5 
Steep        4 

Slope 

Flat          4 
Moderate  2 
Steep        1 

Slope 

Flat         2 
Moderate  0 
Steep        0 

6.8  Clear  cuts,  new  growth  less  than 
3  feet  tall 

10 

Slope 

Flat          7 
Moderate  5 
Steep        4 

Slope 

Flat          5 
Moderate  3 
Steep        2 

Slope 

Flat         3 
Moderate  1 

Steep        0 

6.9  Feedlot  or  concentrated 
livestock  watering 

10 

3 

2 

0 

6.10  Residential  Development 

10 

9 

6 

0 

e.llHuman  constructed  dams  or 

dikes: 

often  indicates  unnatural  wetlands 

Not  Present 
10 

Present 
7 

None  Present 

1-5% 

6-25% 

>25% 

6.12  Human-  induced  saline  seeps 
were  observed 

10 

7 

4 

0 

6.13  Industrial  or  Commercial 
Activities 

10 

7 

4 

0 

6.14  Oil  and  Gas  Development 

10 

7 

4 

0 

6.15  Were  any  of  these  stressors  observed  within  100-  500m  from  the  Wetland?  (Please  circle) 

Row  Crops                                                    Oil  and  Gas  Development                                         Recreational  Activities  (e.g.  campground,  fishing  access  point,  etc.) 
Human-  induced  saline  seeps                            Hayfiek)                                                                Feedlot/concentrated  livestocl<  watering 
Industrial  or  commercial  Activities                        Roads/  Railroad  Grades                                            Clear  cuts  (new  growth  less  than  3  feet  tall) 
Residential  Development                                  Dams  or  Dikes  upstream  (Riverine  Sites) 

Distance  of  road  from  wetland 

>  100  meters 

51-100  meters 

11-50  meters 

<=10  meters 

6.16  2-track  dirt  road  Up  Stope 

10 

6 

4 

2 

6.17  Other  2-track  dirt  road 

10 

8 

6 

4 

6.18Dirtand  gravel  roads,  railroad  grades  Up  Slope 

10 

4 

2 

1 

6.19AII  other  dirt  and  gravel  roads,  railroad  grades 

10 

6 

4 

2 

6.20Paved  Roads  Up  Slope 

10 

2 

1 

0 

6.210ther  Paved  Roads 

10 

4 

2 

1 

Buffer  Condition  Index 

Sum  the  four  lowest  scores  circled  and  divide  by  the  total  nassible  for  the 

Assessment  area  (40).              +         +         +         = 

/40  = 

7.0  ReStOrability  circle  the  appropriate  category  and  sub-category  and  describe  how  the  wetland  is  trending  (when  appropriate) 

7.1  How 

Category  A: 

Category  B: 

Category  C 

Category  D: 

easily  can 

No  observed  impacts; 

Some  slight  impacts  that 

More  significant  impacts  or  disturbances 

Serious  impacts  and  stressors 

the  wetland 

Wetland  does  not  need 

can  be  fixed  or  restored 

within  the  buffer  area  that  can  be  removed. 

are  not  economically  feasible  to 

be  restored? 

to  be  restored. 

with  minimal  expense 

(such  as  a  change  in  land  use  practices: 

remove/restore,  (e.g.,  highway  or 

and  effort  (e.g.  adding 

e.g.  crop  land  changed  to  pasture,  cattle 

fixed  permanent  infrastructure) 

fencing). 

tank  or  abundant  noxious  weeds) 
Restoration  would  require  some  expense 
and  effort. 

7.2  Wetland 

Sub-Category  1: 

Sub-Category  2: 

Sub  Category  3: 

Sub-Category  4: 

Trend 

Wetland  condition  is 

Wetland  condition 

Wetland  condition  is  trending  dovmward. 

Wetland  condition  trend  can  not 

towards 

trending  upward. 

appears  to  be  stable. 

be  determined 

natural 

restoration 

Comments: 
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7.3  Rank  stressor  -  Choose  from  the  ist  and  rank  al  starting  with  1  (hi^est) 

Grazing 

Mining 
Row  Crops 
Rcad/Railroad(s) 
Dam/Dike/Weir 
Extensive  Noxious  Weeds 

Point  Source  Contamination 
Residential  Development 
Human  Recreation 
Industnal  Development 
Forestry/Clear  cutting 

Oil/Gas  Development 
Dredging/Filling 
Feedlot/Cattle  Watering 
De-Watenng 
Hay  Meadow 

Summary  of  Rating 

Hydrooprimnmhir  Dnnriitinn  Inripx          

Veoetatinn  Cnnriitinn  Inriex           

Water  Onalitv  Cnnrlitinn  Index                  

Buffer  nnnrtitinn/.Strf><;<;or  Srnre 

Wetland  Impact  Score  Calculation: 

If  there  is  surface  water  multiply  the  hydrogeomorphic  condition  index  by  0.4;  the  vegetation  condition  index  by  0.4;  the  water  quality 
condition  index  by  0.2. 

If  there  is  no  surface  water  multiply  the  hydrogeomorphic  condition  index  by  0.5;  the  vegetation  condition  index  by  0.5. 


Wetland  Impact  Score. 


Overall  Score  calculations: 


If  there  is  surface  water  multiply  the  hydrogeomorphic  condition  index  by  0.3;  the  vegetation  condition  index  by  0.3;  the  water  condition 
index  by  0.2;  and  the  buffer  condition/  Stressor  index  by  0.2.  Sum  the  indexes  to  detenmine  the  overall  condition  index  score. 

If  there  is  no  surface  water  multiply  the  hydrogeomorphic  condition  index  by  0.4;  tiie  vegetation  condition  index  by  0.4;  tiie  buffer  condition/ 
Sti^essor  index  by  0.2;  Sum  the  indexes  to  determine  ttie  overall  condition  index  score. 


Overall  Score. 


'  This  score  is  not  an  indication  of  wetland  impairment  status.  This  fonn  is  used  to  record  observations  only.  The  form 
can  be  submitted  to  Department  of  Environmental  Quality  for  professional  review  to  assist  in  evaluating  wetland 
condition. 


Overall  condition  index  >0.9-1 .0;  Excellent  Condition 


Overall  condition  index  >0.5-0.7:  Fair  condition 


Overall  condition  index  >0.7-0.9:  Good  Condition 


Overall  condition  index  0.0-0.5;  Poor  Condition 


Appendix  A  -  9 


Appendix  B.  Development  of  the  Montana  Landscape 

Integrity  Model 


The  Montana  Landscape  Integrity  Model  is  designed  to  identify  areas  of  the  state  where  human 
activities  have  changed  the  landscape  in  ways  that  are  broadly  assumed  to  affect  wetland  condi- 
tion by  altering  chemical,  hydrologic,  physical,  or  biotic  properties.    It  uses  four  categories  of 
change:  roads,  land  cover,  hydrology,  and  resource  extraction.  Each  category  has  several  com- 
ponents. For  example,  roads  are  divided  into  four-wheel  drive  roads,  local  roads  and  highways; 
land  cover  is  broken  into  urban  surfaces,  agricultural  lands,  and  timber  harvest  areas',  and  so  on. 

In  ArcGIS.  we  created  an  inverse  distance  raster  layer  for  each  landscape  component,  using  a 
Euclidean  distance  function.  This  resulted  in  a  statewide  data  file  for  each  component,  in  which 
each  30  meter  by  30  meter  pixel  (900  square  meters)  had  a  value  that  represents  its  distance,  in 
meters,  from  the  change  component.  We  then  reclassified  the  data  into  interval  "scores."  based 
on  buffer  distances.  Pixels  within  100  meters  of  a  highway  were  given  a  score  of  5,  for  example, 
while  pixels  more  than  500  meters  away  were  given  a  score  of  1 .  For  large-scale  impacts,  like 
developed  urban  areas,  the  size  of  each  buffer  was  extended;  to  receive  a  score  of  1,  a  pixel  had 
to  be  at  least  2000  meters  from  any  land  cover  classified  as  "developed,  medium  intensity."- 
The  assumptions  underlying  each  score  assignment  were  based  on  the  earlier  data  analysis,  the 
literature,  best  professional  judgment,  or  a  combination  of  all  three.  For  example,  while  we 
assumed  that  the  impacts  of  highways  would  extend  as  far  as  500  meters,  we  surmised  that  the 
impacts  of  four-wheel  drive  roads  would  occur  at  less  of  a  distance,  so  that  data  layer  was  reclas- 
sified into  only  three  levels  (based  on  distances  of  0-100  meters,  100-200  meters,  and  greater 
than  200  meters).  Table  B-1  shows  each  category,  its  components,  and  the  scores  assigned  to 
each  buffer  distance. 

Once  buffer  distances  were  reclassified,  individual  components  were  weighted  and  summed  into 
overall  category  rasters  (e.g.  weighted  roads,  weighted  landcover,  etc)  using  the  weighted  sum 
function.  These  four  category  layers  were  then  weighted  again,  and  combined  into  a  final  Land- 
scape Integrity  Model.  The  road  and  land  cover  categories  were  each  given  35%  of  the  final  total, 
the  hydrology  layer  was  given  20%,  and  the  resource  extraction  category  was  given  a  weight  of 
10%.  Final  values  for  each  pixel  in  the  Landscape  Integrity  Model  ranged  from  100  to  745.  We 
did  not  reclassify  those  values  further  for  this  study,  but  it  could  easily  be  done  to  give  it  a  more 
familiar  scale  (e.g.  1  to  10).  In  this  case,  we  would  also  reverse  the  values  so  that  high  scores 
represented  areas  of  greater  integrity  and  low  scores  represented  areas  of  higher  integrity. 

As  is  always  the  case,  the  model  is  only  as  good  as  the  data  that  goes  into  it.  None  of  the  data 
layers  we  used  in  the  model  is  current  and  complete.  The  roads  layer  captures  section  roads  in 
eastern  Montana  which  are  rarely  if  ever  used  while  omitting  unofficial  but  heavily  used  two- 
tracks,  and  all  the  roads  that  have  been  added  since  the  last  census.  The  land  cover  layer,  while 


'  Timber  harvest  areas  were  identified  from  the  Montana  Cadastral  Layer.  We  selected  lands  owned  by  the  large,  active  timber 
companies:  Plum  Creek,  Stimson.and  Stolze.  and  turned  these  into  a  raster  layer.  Although  timber  harvest  certainly  occurs  on 
smaller  holdings,  we  had  no  way  of  identifying  those  parcels.  Although  the  cadastral  layer  shows  lands  taxed  as  timber  land,  it 
does  not  indicate  what  land  is  being  actively  cut.  When  the  Pacific  Northwest  ReGAP  maps  are  complete,  it  should  be  possible  to 
refine  a  timber  harvest  area 

-  We  did  not  use  the  NLCD  classes  "Developed,  low  intensity"  or  "'Developed,  high  intensity"  in  the  final  iteration.  Low  inten- 
sity development  had  no  significant  relationship  to  wetland  condition  in  our  initial  analysis.  The  NLCD  class  "Developed,  high 
intensity"  is  primarily  a  road  layer  in  Montana,  and  was  therefore  redundant. 
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drawn  from  the  2001  NLCD,  is  based  on  satellite  imagery  from  the  1990s.  The  water  rights 
point  of  use  data  contains  both  georectified  data  with  specific  x,y  coordinates  and  points  dropped 
into  the  center  of  a  public  land  survey  section  because  no  more  precise  data  was  available.  The 
energy  layer  was  current  through  the  summer  of  2008,  but  is  probably  out  of  date  now.  The 
abandoned  mines  data  was  a  point  layer,  and  did  not  necessarily  reflect  the  extent  of  mining 
impacts.  And  some  layers  were  not  available  at  all.    The  absence  of  statewide  NWl  coverage 
meant  that  we  could  not  build  a  distance  layer  from  altered  wetlands,  even  though  the  density  of 
these  wetlands  within  a  given  subwatershed  is  probably  a  very  good  indicator  of  overall  wetland 
condition. 


Table  B-l.  Scoring  and  weighting  of  landscape-level  impacts. 

Buffer  distance 


Category 

(meters) 

Score 

Weight 

Roads 

35% 

4-wheel  drive  (15%) 

0-100 

3 

100.01-200 

2 

>200.01 

1 

Local  roads,  city  streets  (35%) 

0-100 

4 

100.01-200 

3 

200.01-300 

2 

>300.01 

1 

Highways  (50%) 

0-100 

5 

100.01-200 

4 

200.01-300 

3 

300-500 

2 

>500.01 

1 

Land  Cover 

35% 

Urban  (40%) 

0-500 

5 

500.01-1000 

4 

1000.01-1500 

3 

1500.01-2000 

2 

>2000.01 

1 

Crop  agriculture  (40%) 

0-200 

5 

200.01-300 

4 

300.01-400 

3 

400.01-500 

2 

>500 

1 

Timber  harvest  (20%) 

0-500 

5 

500.01-1000 

4 

1000.01-1500 

3 

1500.01-2000 

2 

>2000.01 

1 

Appendix  B-  2 

Category 

Buffer  distance 
(meters) 

Score 

Weight 

Hydrology 

20% 

Artificial  flow  (25%) 

0-100 

3 

100.01-200 

2 

>200.01 

1 

Water  right  point  of  use  (50%) 

0-100 

3 

100.01-200 

2 

>200.01 

1 

Section  404  permit  (25%) 

0-100 

3 

100.01-200 

2 

>200.01 

1 

Resource  extraction 

10% 

Abandoned  mines  (50%) 

0-100 

5 

100.01-200 

4 

200.01-300 

3 

300.01-500 

2 

>500.01 

1 

Oil  or  gas  extraction  (50%) 

0-100 

5 

100.01-200 

4 

200.01-300 

3 

300.01-500 

2 

>500.01 

1 
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Appendix  C.  Ecological  Integrity  Assessment  Form 


SITE  INFORMATION 


SITE  NAME, 
SITE  ID 


ASSESSMENT  AREA  SIZE  IN  M^ 

OWNERSHIP 

HUC4 

HUC5 


DATE  OF  VISIT 
ASSESSED  BY 


PROJECT/PURPOSE 


ELEVATION 

GPSWAYPOINT_ 
Datum_ 

Lat 

Long_ 


Stream  order,  if  rivenne_ 
Fish  sampled'' 


Macroinvertebrates  sampled'. 
Sample  ID,  if  yes 


(Use  decimal  degrees) 


General  site  description,  including  surrounding  uplands 


Directions  to  site; 


Soil  drainage:  Well-drained  Moderately  well-drained 

Total  wetland  area  covered  by  standing  water:  0' 


Poorly  drained  Very  poorly  drained 

'lto25%       26-50%  51-75%         76-100% 


PHOTOS: 


Direction 


Description 


N 


W 


CLASSIFICATION 


ECOLOGICAL  SYSTEM, 
CONFIDENCE  LEVEL: 


Very  High, 


.High, 


Medium 


LOM 


DOMINANT  ASSOCIATION(S):. 


HGM  Wetland  Type:  (Circle  one) 


Riverine 

Upper  Perennial 
Lower  Perennial 
ntermittent 
Ephemeral 


Depressional 

Open 

Closed 

Isolated 


Lacustrine  Fringe 


CONFIDENCE  LEVEL: . 

Comments: 


Very  High, 


-High, 


Medium 


Slope 

Open  Spring 
Rivenne  Spring 
Fen 
Hanging  valley 

Low 


Mineral  Flat 

Wet  Meadow 


COWARDIN  TYPE(S): 
System      Subsystem 


Class 


Water  regime 


Modifier 
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Site  Name 

LEVE 

L  II  ASSESSMENT-Marshes,  wet  meadows,  potholes 

Site  ID 

METRIC 

EXCELLENT(A) 

GOOD(B) 

FAIR  (C) 

POOR(D) 

SCORE 

LANDSCAPE  CONTEXT 

Connectivity 

Non-riverine 

90-100%  natural  habitat 
within  500  m  of  wetland 
penmeter 

60-90%  natural  habitat 
within  500  m  of  wetland 
perimeter 

10-60%  natural  habitat 
within  500  m  of  wetland 
penmeter 

<10%  natural  habitat  within 
500  m  of  wetland  perimeter 

Riverine 

90-100%  natural  habitat 
within  500  m  on  either 
side  and  500  m 
upstream  and 
downstream 

60-90%  natural  habitat 
within  500  m  on  either 
side  and  500  m  upstream 
and  downstream 

1 0-60%  natural  habitat 
within  500  m  on  either 
side  and  500  m  upstream 
and  downstream 

<  1 0%  natural  habitat  within 
500  m  on  either  side  and  500 
m  upstream  and  downstream 

Buffer 

Length 

Buffer  is  >75%  of 
wetland  penmeter 

Buffer  is  50-75%  of 
wetland  penmeter 

Buffer  is  25-50%  of 
wetland  penmeter 

Buffer  is  <25%  of  wetland 
penmeter 

Width 

Average  buffer  width  is 
>  200  m,  adjusted  for 
slope 

Average  buffer  width  1 00- 
200  m,  adjusted  for  slope 

Average  buffer  width  is  50 
100  m,  adjusted  for  slope 

Average  buffer  width  is  <50 
m.  adjusted  for  slope 

Condition 

Buffer  IS  >95%  native 
vegetation  with  intact 
soils  and  little  or  no 
trash  or  refuse 

Buffer  is  >75-95%  native 
vegetation  wnth  intact  or 
slightly  distrubed  soils, 
and  minor  evidence  of 
human  visitation  or 
recreation 

Buffer  is  >25-75%  native 
vegetation  with  slightly  to 
moderately  distrubed 
soils,  and  moderate 
human  visitation  or 
recreation 

Buffer  is  <25%  native 
vegetation  with  severely 
disturbed  soils,  and 
substantial  human  visitation  or 
recreation 

SIZE 

Relative  Patch  Size 

Wetland  is  >95%  of 
onginal  size 

Wetland  is  80-95%  of 
original  size 

Wetland  is  50-80%  of 
original  size 

Wetland  is  <50%  of  original 
size 

Absolute  Patch  Size 

Wetland  is  very  large 
compared  to  others  of 
its  type  (eg.  top  10%) 

Wetland  is  large 
compared  to  others  of  its 
type  (e.g.,  top  10-30%) 

Wetland  is  average 
compared  to  others  of  its 
type  (e.g.,  30-70%) 

Wetland  is  too  small  to 
sustain  full  function  and 
diversity 

VEGETATION  STRUCTURE  (BIOTA) 

Structure 

Vegetation  at  or  near  reference  standard 
condition  in  structural  proportions 

Vegetation  moderately 
altered  from  reference 
standard  condition  in 
structural  proportions 

Vegetation  greatly  altered 
from  reference  standard 
condition  in  structural 
proportions 

Composition 

Vegetation  at  or  near  reference  standard 
condition  In  species  present  and  their 
proportions.  Regeneration  good.  Full  suite  of 
diagnostic  species  present. 

Vegetation  differs  from 
reference  standard 
condition  but  still  largely 
native  Tolerant  or  weedy 
natives  may  be  present 
Many  indicators  absent 

Vegetation  severely  altered 
from  reference  standard. 
Some  strata  absent  or 
dominated  by  weedy  speaes 
Most  indicator  species  absent 

Relative  Cover  of 
Native  Plant  Species 

>99%  relative  cover  of 
native  plants 

95-99%  relative  cover  of 
native  plants 

80-94%  relative  cover  of 
native  plants 

50-79%  relative  cover  of 
native  plants 

Invasive  exotic  species 

No  key  invasive  exotic 
plants  present 

<3%  invasive  exotic 
plants  present 

3-5%  invasive  exotic 
plants  present 

>5%  invasive  exotic  plants 
present 

Organic  Matter 
Accumulation 

Site  has  moderate  amount  of  fine  organic 
matter.  Nev\^  materials  more  prevalent  than 
old  materials  Layers  in  pools  or  topographic 
lows  are  thin. 

Site  is  characterized  by 
small  amounts  of  coarse 
organic  debris,  with  little 
organic  matter 
recuntment,  OR  debns  is 
somewhat  excessive 

Site  has  little  coarse  debris 
and  only  scant  fine  debris  OR 
debris  is  excessive. 

Patch  Types  (See 
below) 

>7  abiotic/biotic  patch 
types  present  in  the 
wetland  {>6  for 
potholes) 

5  to  7  abiotic/biotic  patch 
types  present  in  the 
wetland  (5  or  6  for 
potholes) 

3  or  4  abiotic/biotic  patch 
types  present  in  the 
wetland 

1  or  2  abiotic/biotic  patch 
types  present 

Patch  Interspersion 

Honzontal  structure 
consists  of  a  very 
complex  array  of  nested 
or  interspersed  irregular 
biotic/abiotic  patches 
with  no  single  dominant 
type 

Horizontal  structure 
consists  of  a  moderately 
complex  an-ay  of  nested 
or  interspersed  in-egular 
biotic/abiotic  patches  with 
no  single  dominant  type 

Horizontal  structure 
consists  of  a  simple  an-ay 
of  nested  or  interspersed 
irregular  biotic/abiotic 
patches  with  no  single 
dominant  type 

Horizontal  structure  consists 
of  one  dominant  patch  type 
with  no  interspersion 
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Site  Name_ 
Site  ID 


LEVEL  II  ASSESSMENT-Marshes,  wet  meadows,  potholes 


METRIC 

EXCELLENT(A) 

GOOD(B) 

FAIR(C) 

POOR(D) 

SCORE 

HYDROLOGY 

Water  Source 

Water  source  is 
precipitation, 
groundwater,  natural 
runoff  OR  system 
naturally  lacks  water 
dunng  growing  season 
No  indication  of  direct 
artifical  water  source  or 
point  source  discharge 

Water  source  is  mostly 
natural,  but  site  receives 
occasional  or  small 
amounts  of  inflow  from 
human  sources  e  g  ,  road 
runoff,  storm  drains, 
irrigation)  No  large  point 
source  discharge  into 
site 

Water  source  is  pnmarily 
runoff,  imgation,  pumped 
water,  impounded  water, 
or  other  artificial 
hydrology  Major  point 
sources  discharging  into 
wetland  may  be  present 

Water  flow  has  been 

substantially  diminished  by 
impoundments,  diversions,  or 
withdrawals  from  wetland  or 
adjacent  areas  OR  the  water 
source  is  so  altered  that 
wetland  vegetation  is  gone 

Hydroperiod 

Hydropenod  is 
charactenzed  by  natural 
penods  of 

filling/inundation  and 
drawing  down 

Filling  or  inundation  is 
greater  and  of  greater  or 
lesser  duration  than 
under  natural  conditions, 
but  the  site  is  subject  to 
natural  drying 

Filling  or  inundation  is 
natural,  but  drawdown 
and  drying  more  rapid, 
OR  filling/inundation  is  of 
lower  than  natural 
magnitude  or  duration, 
but  site  IS  subject  to 
natural  drying 

Filling  or  inundation  and 
drawdown/drying  both  deviate 
from  natural  regimes 

Hydrologic 
Connectivity 

Rising  water  in  the  in 
the  site  has  unrestncted 
access  to  adjacent 
upland,  without  levees, 
excessively  high  banks, 
artificial  bamers,  or 
other  obstnjctions  to 
lateral  movement  of 
flood  flows 

Rising  water  has  partially 
restncted  (<50%)  access 
to  upland  due  to 
unnatural  features  OR 
flood  drainage  back  into 
wetland  is  incomplete  due 
to  impoundments  or 
bamers 

Rising  water  has 
significantly  restncted  (50- 
90%)  access  to  upland 
due  to  unnatural  features 

All  water  stages  in  the  wetland 
are  contained  by  artifical 
banks,  levees,  walls,  or  bemns 
or  >90%  of  wetland  has 
bamers  to  drainage  There  is 
essentially  no  hydrologic 
connection  to  uplands 

PHYSIOCHEMICAL 

Soil  Surface  Integrity 

Bare  soil  areas  are 
limited  to  naturally 
caused  disturbances 
such  as  flood  deposition 
or  game  trails 

Bare  soil  due  to  human 
impacts  IS  present  but 
minimal  Water  is  not 
ponding  or  channelled 

Unnatural  areas  of  bare 
soil  are  common  Ponding 
or  channeling  may  be 
present  in  shallow 
disturbances 

Unnatural  areas  of  bare  soil 
are  extensive  and  ponding  or 
channeling  is  likely  Surface 
disturbances  are  deep  and 
widespread 

Water  Quality 

Water  is  clear  with  no 
sheen,  scum,  or  hint  of 
green  Plants  that 
respond  to  enrichment 
are  minimally  present  or 
absent. 

Water  has  a  minimal 
greenish  tint,  cloudiness, 
or  sheen  Plants  that 
respond  to  ennchment 
are  present  but  not 
dominant 

Water  has  a  moderate 
greenish  tint,  sheen,  or 
turbidity  wnth  common 
algae  Plants  that 
respond  to  ennchment 
are  common 

Water  has  a  strong  greenish 
tint,  sheen,  or  turbidity 
Surface  algal  mats  or  other 
vegetation  block  light  to  the 
bottom 

Patch  types: 


Freshwater  marsh 

Open  water-stream 

Oxbow/backwater 

Secondary  channel 

Deep  emergent  plants 

Shallow  emergent  plants 

Beaver  dam 

Trees 

Shrubs 

Spnngs/seeps 

Submerged/floating  veg 

Transitional  meadow 


Pothole 

Open  water 
Shallow  emergent 
Saline  meadow 
Hummocks  or  mounds 
Submerged  or  floating 
Transitional  meadow 
Tall  emergent 
Shrubs 


Wet  Meadow 

Open  water-stream 

Oxbow/backwater 

Secondary  channel 

Deep  emergent  plants 

Shallow  emergent  plants 

Beaver  dam 

Trees 

Shrubs 

Springs/seeps 

Submerged/floating  veg 

Transitional  meadow 


Playa 

Open  water 

Mud  flat 

Salt  flat 

Deep  emergent  plants 

Shallow  emergent  plants 

Saline  meadows 

Greasewood 

Seeps  &  spnngs 

Hummocks  or  mounds 

Submerged  or  floating  vegetation 

Contributing  stream 
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Site  Name_ 
Site  ID 


STRESSORS                                                                                                                                   1 

Land  use  within  SOOm  of  wetland  edge 

Percent  land  use 

Urban  residential 

Industrial/commercial 

Military/airport 

Dryland  farming 

Crop  agriculture 

Orchards/nurseries 

Logging  operation/timber  removal 

Feedlot 

Dairy 

Enclosed  livestock  grazing 

Open  range  grazing 

Sports  field  or  park 

Active  recreation  (OHV,  mountain  biking,  shooting) 

Resource  extraction 

Recent  fire  {<5  years) 

Boating  (motorized) 

Transportation  with  SOOm  of  wetland  edge 

Distance  from  edge 

Lightly  travelled  road 

Moderately  travelled  road 

Heavily  travelled  road 

Pedestrian  trail 

Horse  trail 

Railroad 

Land  use  within  site 

%  of  site 

Mowing 

Livestock  grazing 

Excessive  herbivory 

Excessive  human  visitation 

Tree  cutting/sapling  removal 

Pesticide  or  herbicide  application 

Recent  fire  (<5  years) 

Recent  flood 

Invasive  animals  or  plants 

Hydrology  within  300m 

Impact  (High/Medium/Low) 

Point  source  discharge 

Non-point  source  discharge 

Flow  diversion  or  unnatural  inflow 

Dams 

Flow  obstructions 

Weirs,  headgates 

Dredged  inlet  or  channel 

Engineered  channel 

Dike/levee 

Groundwater  pumping 

Ditches 

Soil  disturbance  witin  SOOm 

Impact  (High/Medium/Low) 

Filling  or  dumping 

Grading/compaction/roadwork 

Plowing  or  discing 

Logging  or  clearing 

Unnatural  areas  of  bare  soil 

Trash  or  refuse 

Pugging,  hummocking,  or  erosion 
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